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INTRODUCTION 

In  a previous  report  [ i]  the  author  presented  a theo- 
retical framework  for  analyzing  the  trajectory  of  rocket- 
satellite  systems,  both  in  the  powered  rocket  ascent  or 
launch  phase  and  in  the  subsequent  orbital  motion  or  coast 
phase.  Techniques  for  solving  the  problem  of  launching 
space  systems  into  proper  earth  orbit  were  elucidated,  and 
subsequently  these  techniques  have  been  tried  on  particular 
cases.  In  the  process,  a numerical  solution  of  the  equa- 
tions of  motion  of  rocket  ascent  has  been  programmed  on  a 
digital  computer  and  found  to  be  useful.  It  is  reported 
here  to  make  it  more  generally  available. 

This  is  obviously  not  the  first  time  a computer  pro- 
gram has  been  developed  to  handle  the  rocket  ascent  pro- 
blem. It  is  possible  to  get  one  of  these  programs  from  an 
outside  source,  and  adapt  it  to  the  specific  tasks  at  hand. 
Very  often,  however,  the  information  accompanying  the  pro- 
gram is  insufficient,  or  the  requirements  of  the  program  are 
not  matched  to  the  data  base  for  the  problem  or  the  com- 
puter at  hand.  It  is  frequently  easier,  faster,  and  more 
reliable  to  develop  a program  ab  initio.  This  was  done, 
and  the  resulting  program  is  presented  and  explained  below. 

It  was  shown  in  Reference  [1]  that  the  rocket  ascent 
problem  could  be  solved  as  if  the  earth  (and  atmosphere) 
were  stationary  in  an  inertial  frame  of  reference,  and 
that  the  corrections  from  earth's  rotation  could  be  added 
in  later.  The  computer  program  developed  here  solves  the 
launch  problem  for  the  stationary  earth  case.  It  will  be 
shown  later  in  a specific  application  how  earth's  rota- 
tion effects  are  adjoined  to  the  computer  program  results. 
For  the  stationary  earth  model  the  simplest  case  of  planar 

rocket  motion  is  shown  in  Fig.  1.  Here  r is  the  radius  vec- 
*► 

Note:  Manuscript  submitted  March  26,  1979. 
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tor  from  the  center  of  the  earth  to  the  center  of  mass  of 
the  rocket,  y is  the  rocket  altitude,  0 is  the  angular 
displacement  from  launch,  which  determines  the  ground  range 
x when  the  radius  of  the  earth  R is  factored  in,  ^ is  the 
flight  path  angle  of  the  center  of  mass  (or  "heading")  of 
the  rocket  with  respect  to  the  local  horizontal,  and  or  is 
known  as  angle  of  attack  of  the  rocket  axis  with  respect 
to  the  flight  direction. 

The  equations  of  motion  for  the  rocket  are  deduced  Cl 3 
straightforwardly  from  the  force  diagram  shown  in  Fig.  2, 
which  also  contains  a brief  description  of  the  indicated 
symbols.  The  forces  indicated  are:  the  gravitational 
force  W,  the  thrust  of  the  escaping  exhaust  gases  T,  which 
is  canted  at  the  angle  fi  to  the  vehicle  centerline,  and  the 
aerodynamic  forces  of  drag  0 and  lift  L due  to  rocket 
motion  through  the  atmosphere.  Drag  and  lift  forces  are 
directed  antiparallel  and  perpendicular  to  the  flight  path 
direction,  respectively.  In  deducing  the  equations  of 
motion,  we  make  the  assumptions  found  in  many  references 
[1-4]: 

(1)  the  aerodynamic  forces  act  through  the  center 
of  pressure, 

(2)  the  force  of  gravity  acts  through  the  center  of 
gravity,  and 

(3)  the  thrust  force  is  applied  through  the  center 
of  combustion. 

The  equations  of  motion  for  the  center  of  mass  motion  are 
determined  [l]  to  be: 

\f  = -(K/r2)  sin  - (D/M)+  (T/fo)cos  («+/3)  (1) 

2 

v $ = - K _y_  cos  + _T  sin  (a+j3)  + __L  (2) 

“r2  r M M 

The  speed  of  the  rocket  center  of  mass  is  here  denoted  by 
v,  and  the  mass  of  the  rocket  is  denoted  by  M.  The  symbol 
K is  the  product  of  the  universal  gravitational  constant 


and  the  mass  of  the  earth,  and  has  the  value 

K = 3.986012X105  kjn/Vsec.2 

The  other  symbols  have  already  been  defined.  The  angle  of 
thrust  £ is  normally  very  small  for  large  launch  vehicles 
(e.g.,  f}/ar«l)  , and  it  is  customary  [2,3]  to  set  fi  to  zero 
in  the  equations  of  motion  before  solving  for  the  launch 
trajectory.  Hence,  the  good  approximation  is  made 

P * 0 (3) 

in  Eq.'s  (1)  and  (2),  thus  avoiding  the  complications  of 
additional  moment  equations  for  the  steering  of  the  rocket. 
The  variation  of  0 is  preprogrammed  in  the  guidance  and 
control  system  of  the  rocket,  which  causes  the  relatively 
much  larger  variations  in  or  involved  in  the  shaping  of  the 
launch  trajectory.  The  approach  to  be  used  in  the  nu- 
merical solution  of  Eq's  (1)  and  (2)  is  in  this  spirit. 

An  ar(t)  profile  will  be  selected  to  cause  the  relatively 
much  larger  variations  in  ^>(t),  and  it  will  be  done  in  such 

a way  as  to  obtain  a correct  value  of  ^ at  rocket  burnout 

and  a smooth  variation  in  $ (t) . The  program  user  typically 
knows  this  burnout  value  of  ^ from  consideration  of  the 
orbit  into  which  he  desires  to  inject  the  payload.  The 
procedure  will  be:  (1)  selection  of  an  a(t)  profile,  (2) 
numerical  solution  of  v(t)  and  ^ ( t ) from  the  equations  of 
motion,  and  then  from  orbit  injection  or  other  criteria, 

(3)  selection  of  a new  a (t)  profile  and  a repetition  of 
step  (2) . This  process  is  iterated  on  the  computer  until 
a "correct"  launch  trajectory  is  obtained. 

Other  approximations  in  the  solutions  of  Eq's  (1)  - 
(3) have  been  experimented  with  during  the  evolution  of 
this  program.  One  approximation  is  the  choice  of  specific, 
analytic  profile  for  *!>(t)  with  free  parameters  adjusted 
so  that  i£(t)  takes  on  specific  end  values.  The  known 
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end  values  are  those  at  rocket  launch  and  at  orbit  injec- 
tion.  The  problem  is  that  all  the  other  values  of  ^(t)  in 
between  are  forced  by  the  analytic  form;  this  will  be  ela- 
borated on  shortly.  The  advantage  of  assuming  an  analytic 
^(t)  profile  is  that  only  Eq.  (1)  needs  to  be  solved  if 
some  other  information  about  or(t)  is  known  - e.g.,  that  a 
is  small  enough  so  that  cos  cr«l  (recall  Eq.  (3)).  The 
criterion  of  small  or  values  is  presumably  desirable  from 
the  standpoint  of  launch  efficiency;  the  rocket  engine  then 
does  its  work  along  the  direction  of  motion.  Ehricke  [3  ] 
uses  this  approach,  but  points  out  that  a could  be  calcu- 
lated from  $(t)  and  the  solution  for  v from  Eq.  (2).  In 
particular,  the  assumption  about  the  smallness  of  a could 
be  checked.  Ehricke  uses  a particular  analytic  form  ip(t) 
for  each  stage  of  a multistage  rocket.  Experience  with 
this  approach,  however,  points  out  certain  flaws  in  it. 

In  the  first  place,  the  end  points  of  ^i(t),  which  nail  down 
its  free  parameters,  are  typically  unknown  ahead  of  time. 
Hence,  there  is  much  time  wasted  experimenting  with  end 
values  (i.e.,  stage  ignition  and  burnout  values  of  ip)  which 
will  keep  the  range  of  a values  small.  A major  problem  is 
that  the  analytic  form  of  ip  ( t)  is  apparently  not  very 
realistic,  because  it  forces  some  rather  wild  variations  of 
a(t).  Neither  is  the  criterion  of  smallness  of  a(t)  neces- 
sarily realistic,  particularly  in  the  higher  rocket  stages 
when  substantial  tilts  may  be  necessary  to  inject  payloads 
at  sufficiently  high  altitudes  to  avoid  aerodynamic  heating 
or  drag  effects.  Ehricke 's  method  was  thus  found  to  be  too 
slow  and  inefficient  in  the  trial  and  error  aspect,  not 
sufficiently  unambiguous  for  inclusion  of  an  automatic 
self-correcting  procedure  based  on  a small  a criterion, 
and  potentially  too  inaccurate  for  our  purposes.  Ruppe[3] 
similarly  employs  analytic  forms  for  sin  ip  and  cos  ip  which 
are  polynomial  expansions  in  t designed  to  give  correct 
launch  and  orbit  injection  end  point  conditions.  With 
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certain  other  approximations  this  enables  one  to  inte- 
grate Eq.  (1)  analytically.  Unfortunately,  the  large  dis- 
agreements in  the  values  inferred  from  the  separate  poly- 
nomial expansions  for  sin  0 and  cos  0 and  the  inaccuracies 
in  the  method  seemed  unacceptable  for  our  purposes. 

Other  approximations,  often  made  in  the  solution  of 
Eq’s  (1)  - (3),  involve  the  evaluation  of  the  force  terms 
in  these  equations.  The  evaluation  of  forces  in  the  pro- 
gram will  be  detailed  in  the  next  section,  but  for  present 
purposes  it  will  suffice  to  just  mention  the  qualitative 
reasoning  behind  the  approximations.  Due  to  ambient  pres- 
sure decrease  with  altitude  it  is  found  that  the  effective 
rocket  thrust  force  T in  Eq's  (1)  and  (2)  increases  from 
its  launch  value  Tgl  at  sea  level  to  its  vacuum  value  Tvac 
at  the  orbit  injection  altitude.  This  increase  of  thrust 
effect  is  at  least  partially  offset  by  the  drag  force  in 
Eq.  (1)  which  acts  to  slow  the  rocket  down.  This  has 
prompted  certain  authors  [3,4]  to  eliminate  the  drag  and 
thrust  increase  effects  from  Eq.  (1) , since  they  tend  to 
counterbalance  each  other,  and  use  a constant  thrust  term 
instead.  This  is  a nice  simplification,  particularly  since 
the  drag  and  lift  coefficients  for  the  rocket,  which  enter 
in  the  evaluation  of  D and  L in  Eq's  (1)  and  (2),  are 
usually  not  known  very  well.  Similarly,  the  force  L is 
often  dropped  from  Eq.  (2)  on  the  grounds  that  it  is  not 

very  important  for  large  rockets.  It  has  been  the  author's 
experience  that  for  l^rge  rockets,  e.g.,  with  liftoff 
weights  of  hundreds  of  thousands  of  lbs.  or  more,  the  thrust 
increase  effect  outweighs  the  drag  effect.  As  a result, 
the  constant  thrust  approximations  give  injection  velo- 
cities smaller  than  the  more  complete  evaluations  by 


as  much  as  2-3%.  This  kind  of  discrepancy  is  serious  for 
the  orbit  determination,  as  shown  in  Fig.  3,  which  gives 
the  dependence  of  orbit  apogee  altitude  on  perigee  injec- 
tion velocity.  A two  percent  variation  injection  speed 
can  make  a difference  of  500  nmi.  in  apogee  altitude.  In 
questions  of  orbit  determination  (or  even  whether  or  not  a 
stable  orbit  has  been  achieved)  this  kind  of  uncertainty 
may  well  be  unacceptable.  As  Eq.  (1)  indicates,  the  thrust 
term  also  varies. when  the  angle  of  attack  a becomes  large, 
which  is  often  the  case  in  the  upper  rocket  stages.  The 
constant  thrust  approximation  was  therefore  abandoned.  The 
decision  was  also  made  to  not  discard  the  lift  term  L in 
Eq.  (2),  since  it  was  found  [3]  to  significantly  enhance 
the  leverage  exerted  by  the  assumed  a(t)  profile  on  the  cal- 
culated j[i(t)  solution  at  ascent  altitudes  where  aerodynamic 
effects  attain  their  largest  values. 

For  the  above  reasons  the  approach  adopted  in  the  com- 
puter program  is  to  assume  an  or ( t ) profile,  solve  the  equa- 
tions of  motion,  and  then  refine  the  original  or(t)  in  an 
iterative  approach.  This  proved  to  be  easy  to  automate, 
and  the  results  did  not  appear  to  be  sensitive  to  the  am- 
biguities in  choosing  the  ar(t)  profile.  With  an  evaluation 
of  thrust  increase  and  aerodynamic  effects,  and  with  a non- 
reliance on  the  constant  thrust  approximation,  the  program 
has  reached  a level  of  approximation  where  it  should  be 
useful  for  most  purposes,  including  orbit  determination. 
More  on  this  later,  but  the  next  task  is  to  supply  details 
about  the  program  and  its  use. 

The  Evaluation  of  Terms  in  the  Equations  of  Motion 

For  a rocket  with  swivel  control  motors  the  total 
thrust  T in  Eq's  (1)  and  (2)  includes  pressure  forces,  and 
is  written  [3-5]  as 


-6- 


T = Mve  + A [pe  - p(y)3  / (4) 

where  M is  the  mass  loss  rate  of  exhaust  gases,  v is 
the  actual  average  axial  speed  of  these  gases  relative  to 
the  rocket,  A is  the  exit  aperture  area,  pQ  is  the  pres- 
sure of  the  exhaust  gases  at  the  exit  aperture,  and  p(y) 
is  the  ambient  pressure  of  the  atmosphere  at  altitude  y. 

This  pressure  varies  from  p(o)  to  zero  during  the  first 
part  of  rocket  flight,  usually  by  the  time  of  first  stage 
separation.  Consequently,  T varies  from  its  sea  level  value 


T , to  its  vacuum  value  T 
si  vac 


Both  T . and  T are  nor- 
sl  vac 


mally  specified  for  the  first  stage  of  the  rocket,  pos- 
sibly also  for  the  second  stage,  £>ut  usually  only  Tv&c  is 
given  for  higher  stages.  Hence,  for  the  first  stage,  pos- 
sibly also  the  second,  the  ratio 

x = T /T  . (5) 

vac  si 

is  known,  and  Eq.  (4)  can  be  rewritten  as 

T-T  P(^.._l  (6) 

L x p (o)  _ 

The  computer  program  evaluates  thrust  from  this  equation, 
which  exhibits  the  increase  of  thrust  with  altitude. 

The  aerodynamic  forces  of  lift  L and  drag  D in  Eq's 
(1)  and  (2)  can  be  evaluated  from  [3] : 


D = [CDo  + CDLa2]  <Pv2/2)S 

L -[(dqydot)  or]  (Pv2/2)S 


where  p is  the  atmospheric  mass  density,  v is  the  rocket 
speed  and  S is  a reference  cross-sectional  area  (e.g., 
that  for  the  lowest  rocket  stage)  to  which  the  aerodynamic 
coefficients  are  referred.  As  seems  plausible,  the  aero- 
dynamic forces  are  seen  to  be  proportional  to  the  dynamic 
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pressure  of  the  atmosphere  (i.e.,  pv2~/2)  and  the  reference 
cross-sectional  area  of  the  rocket  vehicle.  The  propor- 
tionality constants  in  square  brackets  in  Eq(s  (7)  and  (8) 
are  the  drag  and  lift  coefficients  which  depend  on  both  the 
angle  of  attack  or  and  the  speed  v.  The  dependence  on  a is 
to  be  expected,  since  the  tilting  of  the  rocket  centerline 
away  from  the  flight  direction  obviously  exposes  more  of 
the  rocket  surface  to  the  pressure  forces  of  the  atmosphere, 
thus  increasing  the  aerodynamic  forces.  The  lift  force 
vanishes  for  or  = o.  The  representation  of  the  a dependence 
of  the  drag  and  lift  coefficients  in  Eq's  (7)  and  (8)  is 
supposed  to  be  valid,  according  to  Ehricke  [3],  for  moder- 
ate angles  of  attack  (e.g.,  a < 10°).  Here  C„,  C__  , and 
dC_/d  a are  constants  which  depend  only  on  the  speed  of  the 
rocket.  Actually,  the  dependence  is  found  to  be  on  the 
mach  number  M,  defined  as 

M = v/c  , (9) 


where  c is  the  speed  of  sound  at  the  altitude  of  the  rocket. 
We  have  used  the  dependences  given  by  Ehricke  (cf.  his  Fig. 
5-8  in  [3])  for  a two-stage  rocket,  a close  analytic  fit 
to  which  has  been  found  to  be: 


0.25  + 0.3773  exp  [ -16 . 6426 (M-l. 14) 2] 

0.6273  - 30.902  |m-1.14|3,5 
0.17913  + 0.5318/M 


(OSM*  0.88889) 
(0.88889<M£1.28) 

(i.28  < m)  (io; 


c 1.822  - &298  cos  [M  tt/1.086] 
DL  j-0.2818  + 3.2727/M 

dCj/da  = { 0.8  + 6/M 


(0  s 1.5) 
(1.5  < M) 


(11) 


(0  C M S |.5) 

(M  > 1.5)  (12) 


The  graphs  of  Eq's  (10)  - (12)  are  shown  in  Fig.  4,  along 
with  a sketch  of  the  type  of  rocket  to  which  these  coef- 


ficients apply.  In  this  case,  the  reference  value  S in 


Eq's  (7)  and  (8)  is  the  cross-sectional  area  of  the  bottom 
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rocket  stage.  The  dynamic  pressure  factor  in  Eq 1 s (7)  and 
(8)  is  seen  to  involve  a competition  between  an  atmos- 
pheric density  p which  is  decreasing  exponentially  with 
altitude  and  a quadratic  dependence  on  velocity  which  is 
increasing  with  altitude.  The  result  of  this  in  calcu- 
lations is  that  aerodynamic  forces  are  negligible  above 
about  60  km,  and  peak  somewhere  «round  10  km. 

For  the  evaluation  of  thrust  and  aerodynamic  forces 
an  atmospheric  model  is  needed,  i.e.,  a specification  of 
pressure  p and  density  p vs  altitude  y in  the  region  in 
which  thrust  variation  and  aerodynamic  forces  play  a role 
(e.g.,  for  0 £ y < 100  km).  The  temperature  is  related  to 
p and  p by  the  ideal  gas  law; 


p = p RT/M  (13) 

where  M is  the  mean  molecular  weight  of  one  mole  of  atmos- 
pheric gas,  and  R is  the  gas  constant  given  by  R = 8.3143 
joules  mole  ^deg  8.3143X  lO^ergs  mole^deg  The 

speed  of  sound  in  Eq.  (9)  is  calculated  from 

c^  = 1.4  p/p  (14) 


where  the  factor  1.4  is  the  ratio  of  specific  heat  at  con- 
stant pressure  to  the  specific  heat  at  constant  volume  for 
diatomic  molecules,  such  as  air.  For  use  in  the  computer 
program  it  would  be  helpful  to  have  an  analytic  expression 
of  the  atmospheric  density  and  pressure  altitude  profiles. 
This  is  developed  in  Appendix  A.  Using  the  approach  de- 
veloped there  and  the  numerical  tabulation  of  Jastrow  and 
Kyle  [6],  the  following  analytic  forms  for  the  atmospheric 
profile  are  obtained; 

p = po(H/Ho)_1/^o 

P = pn  [1-2  (y-y_)/(R+y  )]-1  (h/h  ) “(1+?o} /po  ' 


r 


where  h = hq  + £o(y-yo)  (y^y^Yf)  (15 

Different  parameters  ( pQ , pQ,  Hq,  fiQ,  yQ)  are  used  for  dif- 
ferent altitude  intervals  [Y^  YfL  and  these  are  given  in 
Table  1.  Use  of  these  expressions  is  found  to  give  a good 
fit  of  the  numerical  tabulation  of  Jastrow  and  Kyle  [6]. 
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Integration  of  the  Equations  of  Motions 


Before  going  on  to  the  actual  integration  of  the  equa- 
tions of  motion,  further  approximations  used  should  be  in- 
dicated. These  have  to  do  with  the  magnitude  of  the  angle 
of  attack  at  during  the  flight  of  the  rocket. 

As  Eq.  (2)  indicates,  while  the  rocket  speed 
is  relatively  low  during  first  stage  motion,  a particular 
value  of  a has  a larger  effect  on  the  turning  rate  p than 
when  the  rocket  speed  is  much  greater  in  second  or  higher 
stage  motion.  Values  of  a in  first  stage  motion  needed  to  a 
cause  a particular  orbit  injection  value  of  ip  (e.g.,  ip=0) 
have  been  found  not  to  exceed  a few  degrees.  On  the  other 
hand,  the  maximum  value  of  a in  higher  stages  has  been  found 
to  be  as  high  as  tens  of  degrees  in  some  cases,  depending  on 
payload  weight  and  the  required  orbit  injection  altitude. 

At  least  while  aerodynamic  effects  are  important,  however, 
it  appears  to  be  a good  approximation  to  replace  sin  a by  a, 
and  the  approximations  of  Eq's  (7)  and  (8)  are  valid.  Hence 
to  a good  approximation,  Eq's  (1)  - (3)  become 

v = -(K/r2)  sin  ip  - (D/M)  + (T/M)  cos  a (16) 

v$>  = -[k  - v 1 cos  ip  + 1 [*T  + SCL  Pv  s 

p-  — J M [ ^ T 


Is  in  or 

(17) 


These  are  the  equations  of  motion  integrated  in  the  computer 
program.  They  are  considered  valid  for  all  rocket  stages. 
Here  D is  given  by  Eq's  (7),  (10),  and  (11),  and  T is  given 
by  Eq.  (6) . 

The  procedure  is  to  employ  a particular  profile  for 
or (t ) in  the  integration  of  Eq's  (16)  and  (17).  It  has  the 
form 


C (t-t  ) [l-  (t~to) 
l td 


a(t)  = 


(to<t<Vtd) 


(18) 


It  is  shown  in  Fig.  5 as  an  inverted  parabola  where  it  is 
non -vanishing,  with  an  extremum  point  at  t = tQ  + t^/2  and 
associated  value  arm  = Ct^/4.  It  is  a convenient  profile  to 
manipulate,  and  the  results  for  orbit  injection  conditions 
have  been  found  to  be  insensitive  to  variation  of  para- 
meters in  the  profile.  For  example,  a particular  orbit  in- 
jection altitude  at  ip  = 0 can  be  obtained  for  various  C 
values  by  adjusting  the  td  value,  and  the  value  of  injec- 
tion velocity  is  insensitive  to  this  procedure  (e.g., 
within  0.2%).  The  assumption  of  a particular  form  in  Eq. 

CL8)  would  be  a major  concern  if  this  were  not  the  case. 

Beside  having  solutions  for  v(t)  and  ^i(t)  in  Eq's  (16) 
and  (17),  it  would  also  be  of  interest  to  find  the  altitude 
and  ground  range  (cf . Fig.  1) , given  as  solutions  of 

9 = v sin  ip  (i9)" 

* = v R cos  ip/( R+y)  (20) 

The  Eq's  (16),  (17),  and  (19),  and  (20)  are  a coupled  set 

of  ordinary  first  order  differential  equations  in  the  time 
variable,  and  are  therefore  solvable  by  the  Runge-Kutta 
method  [7].  In  explaining  this  method  it  will  be  convenient 
to  rewrite  these  equations  as 

x(a)  = f(«)  (t,  Cx(^)])  (®, 0=1, 2,3,4)  , (21) 

where  the  superscript  a in  this  case  runs  over  the  four 
variables  v,  ip,  y,  and  x,  which  are  respectively  repre- 
sented by  x(1>,  x and  x^.  Eq.  (21)  states 

that  the  time  derivative  of  one  of  these  variables  is  a 
particular  function  of  the  time  (e.g.,  through  the  given 
explicit  time  dependences  of  or  ( t ) in  Eq.  (18)  and  mass  of 
the  rocket  M(t))  and  of  the  four  variables  {x^3  at  that 
time.  Actually,  inspection  of  Eq's  (16),  (17),  (19),  and 
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— 


* 


on  x 


(a) 


( 4 ^ = x^ere.  In  the 


(20)  reveals  no  dependence  of  f'”'  on  x ’ = x 
Runge-Kutta  numerical  method,  the  time  scale  is  broken  up 
into  mesh  points  separated  by  a fixed  interval  width  h. 

The  integration  is  consecutively  forward  on  the  time  scale, 
with  the  variable  at  the  mesh  point  n determined  from  the 
known  values  at  the  mesh  point  n-1,  as  follows: 


(a) 

xn 

(a) 

" xn-l 

+ 1 (doc.  (a)+  2 Ax 
6 1 

where 

= hf  ^ 

(t  t x 

'n-1 * 1 xn-l  J 

= hf(0,) 

(tn-l  + h/2'  txi 

= h t{a) 

(tn_i  + h/2 , [ Xj 

*C4W 

= hf'"* 

(tn-l+h'tx„-^ 

24x^0r)  +Ax4(0f))  (22) 


Ax/^/2}) 


(23) 


The  procedure,  according  to  these  equations,  is  to  first 
evaluate  the  four  values  of  £x,  ^ , which  then  enables  the 
evaluation  of  the  four  values  of  Ac* [ ' . Next  comes  the 
evaluation  of  and  then  . Finally,  xn^  is 

evaluated  from  Eq.  (22) . This  is  also  the  procedure 

for  any  number  of  variables  («,  f3,  = 1,2,  ),  not  just 

the  four  of  interest  here,  and  the  method  is  easily  pro- 
grammable . 

The  computer  program  for  the  rocket  trajectory  con- 
sists of  a main  program  and  ten  subroutines.  The  general 
approach  of  the  coordinating  main  program  RALLY  6 is  to 
take  a trial  ar(t)  function  from  Eq.  (18),  and  then  perform 
the  numerical  integration  of  Eq's  (16),  (17),  (19),  and 
(20)  by  the  Runge-Kutta  method  outlined  in  Eq's  (21)-(23). 
Based  on  the  final  rocket  burnout  values  of  ^(t)  and  y(t) 

- i.e.,  the  orbit  injection  values  of  these  variables  - 
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the  tilt  parameters  in  Eq.  (18)  are  varied,  and  the  pro- 
cess is  repeated  iteratively  until  the  desired  orbit  in- 
jection conditions  are  obtained.  Instead  of  varying  tilt 
parameters  the  program  also  provides  for  the  variation  of 
coast  times  after  powered  stage  burnouts,  in  order  to  obtain 
orbit  injection  conditions  . The  program  is  interactive  with  the 
user  in  that  corrections  are  made  by  him  during  the  run, 
but  at  some  point  the  user  can  choose  an  option  which 
makes  the  program  self-correcting,  i.e.,  the  program  iter- 
ates on  the  corrections  by  a Newton -Raphson  method  [7]  un- 
til convergence  to  a specified  (by  the  user)  stage  burnout 
heading  is  obtained  within  a specif ied  tolerance. 

The  input  and  output  options  are  very  flexible,  so  that  the 
program  will  handle  a wide  variety  of  rocket  problems  with 
ease.  The  Subroutines  are:  (1)  MASS  which  evaluates  the 
mass  of  the  rocket  in  Eq's  (16)  and  (17)  from  M=Mq -Mt, where 
the  initial  stage  mass  Mq  and  mass  ejection  rate  M are 
known  parameters  for  the  rocket  stage,  and  t is  the  burn 
duration  in  that  stage;  (2)  TILT  which  evaluates  angle  of 
attack  from  Eq.  (18) ; (3)  DEN  which  evaluates  density  and 

pressure  according  to  Eq.  (15)  and  Table  1 and  sound  speed 
according  to  Eq.  (14) ; (4)  PSINCR  which  evaluates  ^ in  Eq. 

(17) , which  is  needed  for  the  Runge-Kutta  integration  in 
the  main  program; (5)  VINCR  which  evalues  v in  Eq.  (16) ; 

(6)  ALT  which  evaluates  y in  Eq.  (19) ; (7)  GRRG  which 
evaluates  x in  Eq.  (20) ; (8)  CNVG  which  computes  the  para- 
meter corrections  to  tilt  amplitude  (C  in  Eq.  (18))  or 
stage  coast  time  in  the  self-correcting  mode  of  operation 
of  the  program,  and  then  tests  for  convergence;  (9)  STORE 
which  either  stores  burnout  values  just  computed  for  the 
present  rocket  stage  to  be  used  in  calculations  for  the 
next  stage,  or  recalls  values  stored  from  calculations  on 
the  preceding  stage  to  be  used  in  a repetition  of  the  cal- 
culations for  the  present  stage;  and  (10)  FORCES  which  com- 
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putes  thrust  and  aerodynamic  forces  in  the  equations  of 
motion  according  to  Eq's  (5) -(12). 

Detailed  Description  of  the  Computer  Program  and  Its  Use 
In  Appendix  B is  the  listing  of  the  program  for  the 
DEC-10  system,  and  in  Appendix  C is  the  almost  identical 
listing  of  the  program  for  the  PDP-11/70  computer  system. 

One  of  the  few  changes  made  arose  from  the  way  the  two 
different  computers  evaluated  cos  # in  Eq's  (17)  and  (20). 
There  was  no  problem  with  the  DEC-10,  but  the  PDP-11/70 
made  small  errors  in  the  evaluation  of  cos  $ near  ip=n/2 
rad.  As  small  as  these  errors  were,  they  were  compounded 
in  the  integration  procedure,  so  that  final  errors  in  the 
burnout  ground  range  and  altitude  amounted  to  a few  km. 
or  more.  This  was  remedied  in  the  PDP-11/70  program 
version  by  replacing  cos  0,  wherever  it  occurred,  by 
sin  (tt/2-$)  . With  respect  to  the  DEC-10  listing  it  is 
seen  that  statements  #12400  and  14900  were  changed  in  the 

PDP-11/70  version  in  this  fashion.  There  is  also  a dif- 
ference in  the  treatment  of  the  vacuum-to-sea  level  thrust 
ratio  which  is  given  a value  VSLR(I)  for  each  rocket  stage 
in  the  DEC-10  version  (cf.  statements  #120,  2010,  and  4255 
in  Appendix  B) . This  ratio  is  given  as  x in  Eq's  (5)  and 
(6) . Very  often  the  effect  of  VSLR  or  x is  negligible,  ex- 
cept for  the  first  stage,  because  rocket  altitude  in  the 
upper  stages  is  high  enough  that  the  atmosphere  is  essen- 
tially a vacuum,  i.e.,  T takes  on  the  value  T in  Eq.  (6)  . 

VflC 

This  is  why  VSLR  is  treated  as  a single  parameter  for  the 
first  stage  in  the  PDP-11/70  version  in  Appendix  C.  The 
only  other  change  was  necessitated  by  the  difference  in 
printing  output  between  the  two  systems.  While  the  output 
is  typed  at  the  DEC-10  terminal  on  a paper  roll  as  TYPE 

statements  are  encountered  in  the  program,  it  is  printed  on 
a CR-tube  display  in  the 
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PDP-11/70  system  (which  periodically  erases  itself)  and 
subsequently  printed  on  paper  at  a line  printer  terminal 
upon  successful  completion  of  the  job.  In  the  PDP-11/70 
version  of  the  program  the  IPR  printing  option  variable  in 
the  DEC-10  version  is  replaced  by  ITY,  and  IPR  becomes  a 
printing  option  variable  for  the  line  printer.  Otherwise, 
the  programs  are  identical.  Henceforth,  the  discussion 
will  be  confined  to  the  DEC-10  version. 

At  this  point  the  statements  of  the  program  listed  in 

Appendix  B will  be  itemized#  and  to  help  the  reader  follow 

the  logic  of  the  program,  a flow  chart  of  the  main  program 

(RALLY  6)  is  included  in  Fig.  6(a),  (b) , and  (c) . The  first 

two  DATA  statements  in  Appendix  B (#300  and  400)  list  the 

yf  and  yQ  values  in  Table  1,  while  the  next  one  lists  the 

values  of  pQ  in  g./cm.3.  The  next  data  statement  (#700) 

lists  the  values  of  -2(R+y0f3in  km-1  (cf.  Eq.  (15)).  Next 

come  the  values  of  (in  #1000),  (in  #2000),  and  P in 
2 ° ° o 

dynes/cm  (in  #1350)  from  Table  1.  This  is  followed  by 

values  for  constants  used  in  the  program;  sequentially, 

3/2 

these  are  K in  km  sec  (cf  Eq's  (1)  and  (2)),  the  radius 
of  the  earth  R in  km,  pi,  and  the  pressure  at  sea  level 
(p(o)  in  Eq.  (6))  in  dynes/cm  . After  zeroing  a few  exe- 
cution option  parameters,  the  user  reads  in  the  number  of 
rocket  stages  NS  up  to  6 (which  can  be  changed  to  a higher 
number,  if  necessary,  by  changing  statement  #120)  and  SMZ 
(I),  the  initial  mass  for  each  rocket  stage  (in  lbs.), 
which  is  defined  as  the  mass  of  the  rocket  at  the  begin- 
ning of  the  stage.  Periods  of  coasting  can  be  included  as 
separate  stages.  Next  SMD  (I)  is  read-in  which  is  the 
rate  of  mass  loss  of  each  stage  due  to  fuel  consumption  (in 
lbs/sec).  This  would  be  set  to  zero  for  a coasting  stage. 
Then  the  vacuum  thrust  for  each  stage  TH(I)  is  read  (in  lbs), 
which  is  the  Tvac°f  Eq.  (6)  . Following  this,  the  program 
calls  for  the  vacuum-to-sea  level  thrust  ratio  VSLR  for 
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each  stage  (i.e.,  x in  Eq.  (5)),  the  duration  time  (in  sec) 
for  each  stage  TBD(I),  the  stage  diameter  DI(I)  (in  meters), 
and  the  number  of  time  intervals  NI  (I)  that  each  rocket 
stage  is  to  be  broken  up  into  in  the  specification  of  the 
time  mesh  for  the  numerical  integration.  In  statements 
#2605  - 2630,  TBS (I)  is  computed,  which  is  the  time  mea- 
sured from  launch  in  sec.  at  which  the  I'th  rocket  stage 
commences.  Next  the  altitude  YI,  time  TI,  speed  VI, 
heading  PS  I,  altitude  interval  ID  (associated  with  Table  1) 
and  angle  of  attack  AL  are  given  values  appropriate  to 
launch  conditions,  and  the  program  is  then  ready  to  enter 
the  numerical  integration  procedure. 

After  testing  a pair  of  parameters,  which  will  be  dis- 
cussed later,  the  program  calls  for  the  specification  of 
the  angle  of  attack  parameters  in  statements  #3460  - 3550. 
Here  CA  TLZ , and  TLD  refer  to  C (in  rad.),tQ  (in  sec),  and 
t^  (in  sec)  in  Eq.  (18),  and  thus  refer  to  the  tilt  ampli- 
tude, the  time  measured  from  launch  at  which  tilt  begins, 
and  the  duration  of  non-vanishing  tilt,  respectively.  Now 
the  user  supplies  the  parameters  which  determine  the  various 
execution  options  in  statement  #3575: 

(1)  NT  specifies  the  number  of  rocket  stages  ♦ ' be 
integrated  as  in  Eq’s  (21)  - (23). 

(2)  IOP  specifies  what  part  of  the  program  that  control 
is  returned  to  after  the  NT  rocket  stages  have 
been  integrated.  IOP  = 1 terminates  the  program, 

IOP=2  returns  control  to  program  statement  3 (PS3 ) 
where  variables  are  given  their  launch  values.  IOP=3 
returns  control  to  PS69  which  redefines  the  time  mesh 
for  the  numerical  integration.  IOP  = 4 returns  con- 
trol to  PS25  which  initiates  the  integration  of 
rocket  stages  other  than  the  first.  I0P=5  returns 
control  to  PS2  which  defines  a new  rocket  problem. 

(3)  ICE  specifies  the  number  of  rocket  stages  which 


have  previously  been  integrated. 

(4)  IPR  specifies  the  extent  of  the  printed  output  from 
the  numerical  integration  of  the  NT  rocket  stages. 

IPR  = 0 causes  printout  of  the  variables  at  all  the 
mesh  points.  IPR  fe  1 causes  printout  of  only  the 
final  or  burnout  values  of  the  variables.  IPR  = 0.1 
signals  the  program  that  corrections  to  tilt  or  coast 
time  parameters  will  be  supplied  by  the  user.  IPR=2 
signals  that  tilt  amplitude  or  coast  time  corrections 
will  be  computed  automatically  by  the  program  in  its 
self -correcting  mode  until  convergence  has  been  ob- 
tained, i.e.,  when  the  burnout  value  of  the  heading 
(0)  takes  on  a specified  value  within  a specified 
tolerance. 

(5)  IBG  specifies  what  the  value  of  IOP  will  be  after 
convergence  has  been  obtained  in  the  self-correcting 
mode  of  the  program. 

(6)  JCV  specifies  what  is  to  be  corrected  in  the  self- 
correcting  mode  of  the  program.  JCV  = 2 specifies 
that  TBD(NS)  is  to  be  adjusted,  i.e.,  the  coast  time 
of  stage  NS.  JCV  * 2 specifies  that  CA  is  to  be  ad- 
justed, i.e.,  the  amplitude  parameter  for  the  angle 
°f  attack  profile  of  a particular  rocket  stage. 

(7)  IMA  specifies  whether  or  not  the  payload  weight  is  to 
be  adjusted.  This  control  is  convenient  for  obtaining 
the  dependence  of  orbit  injection  conditions  on  pay- 
load  weight. IMA  - 1 specifies  that  a payload  weight 
increment  is  to  be  read-in  after  the  program  has  ob- 
tained covergence  in  its  self-correcting  mode  of 
operation. 

(8)  JNS  specifies  whether  or  not  a new  NS  is  to  be  read 
in  for  the  integration  procedure.  JNS  = 0 says  not, 
but  JNS  ^ 0 specifies  that  a new  NS  will  be  read  in. 

A temporary  NS  value  is  used  for  the  self-correcting 
mode  of  operation  on  parameters  of  an  intermediate 
rocket  stage  when  the  NT  stages  being  integrated  do 
not  include  the  final  rocket  stage. 

If  IPR  = 2 the  program  enters  its  self-correcting 

mode  of  operation  in  the  integration  of  NT  rocket  stages. 

The  first  part  of  this  is  the  specification  of  the  desired 
final  result  PSSS  for  the  heading  in  the  integration  and  the 
convergence  within  EPS  of  the  parameter  that  is  being  cor- 
rected. These  values  are  read  in  statements  #3882  and  3885. 
The  program  then  carries  out  the  numerical  integration  of 
the  Eq's  of  motion  (16),  (17),  (19),  and  (20)  as  indicated 


in  Eq's  (21)  - (23),  and  with  the  help  of  the  subroutines 
discussed  at  the  end  of  the  preceding  section.  The  reader 
will  have  little  difficulty  in  following  this.  While  some 
of  the  parameters  were  read  in  CGS  or  English  units,  every- 
thing is  converted  to  MKS  units  in  the  numerical  integration 
with  the  exception  that  distance  units  are  kn.  in  the  final 
results . 

Beginning  with  statement  #8010  the  program  enters  a 
series  of  tests  which  are  designed  to  correctly  set  the 
parameter  ICD  for  the  subroutine  STORE,  which  is  called  in 
statement  #8050.  If  ICD  = 1,  STORE  will  store  the  tra- 
jectory values  just  computed,  to  be  used  as  initial  values 
for  integration  of  the  next  stage.  If  ICD  = 2,  STORE  will 
recall  the  initial  values  for  the  integration  just  computed, 
so  that  this  integration  can  be  repeated  to  get  a heading 
which  is  closer  to  that  desired.  The  rest  of  the  main  pro- 
gram, which  is  indicated  by  the  flow  chart,  is  mostly  con- 
cerned with  obtaining  convergence  in  the  self -correcting 
mode  of  operation.  Convergence  is  on  the  tilt  amplitude 
CA  or  on  the  coast  time  of  the  stage  NS  (i.e.,  TBD(NS)),  as 

determined  by  JCV.  As  previously  stated,  the  sub-routine 
CNVG  computes  the  corrections  to  these  parameters  in  the 
self -correcting  mode  (see  statement  #8655).  The  parameter 
ICVG  signals  whether  or  not  convergence  has  been  obtained. 
ICVG  = 0 means  it  has  not;  ICVG  = 1,2  means  it  has. 

Perhaps  the  best  way  to  understand  how  the  program 
works  and  how  to  use  it  is  to  trace  through  the  flow  chart 
in  Fig.  6 with  data  on  sample  problems.  For  this  purpose, 
consideration  is  directed  to  a two-stage  rocket  with  four 
stages  of  motion.  The  author  would  like  to  make  the  dis- 
tinction between  rocket  stages  and  stages  of  motion.  The 
program  is  solely  concerned  with  the  latter  usage,  and  the 
subsequent  discussion  here  is  too.  The  first  stage  of 
motion  is  powered,  the  second  stage  is  defined  as  a 
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coasting  stage,  the  third  stage  is  powered,  and  the  fourth 
stage  is  a coasting  stage  for  the  payload.  It  is  assumed 
that  all  the  rocket  data  have  been  read  in,  so  that  the  point 
in  the  program  has  been  reached  where  the  angle  of  attack 
dataare  to  be  read-in,  ije.,  at  statement  #3500  in  the 
program.  The  solutionsof  four  problems  for  this  rocket  are 
of  interest,  and  these  are  given  below. 

Problem  1 ; The  right  amount  of  tilt  must  be  put  in  the 
first  stage  so  that  the  rocket  will  be  injected  into  orbit 
at  zero  heading  (i.e.,  = 0 as  the  final  solution  of  the 

equations  of  motion) . Hence,  the  second  through  fourth 
stages  fly  in  a gravity-turn  (o=o.inFig,  1)  trajectory.  The 
user  must  supply  the  data  shown  in  Table  2 , which  also  shows 
the  statement  number  that  demands  the  data, along  with  occa- 
sional comments.  The  reader  should  trace  the  path  of  logic 
through  the  flow  chart  which  is  dictated  by  the  data.  The 
initial  reading  of  the  tilt  or  angle  of  attack  data  will 
result  in  a heading  which  is  non-zero  in  the  printed  out- 
put. This  suggests  a new  CA  value  which  is  read-in  (see 
Table  2) , and  the  change  of  IPR  from  1 to  2 in  the  next 
data  input  will  cause  the  program  to  go  into  its  self- 
correcting  mode  cf  operation.  From  then  on  the  program 
iterates  on  CA  until  corrections  to  CA  become  less  than 
0.000002.  This  will  typically  take  only  a few  seconds  of 
actual  runtime  on  the  DEC-10  system.  Then  the  complete 
launch  trajectory  is  printed  out  at  every  time  mesh  point 
as  the  converged  solution.  Note  that  the  user  supplies 
two  trial  values  of  CA  before  the  correction  procedure  is 
made  automatic.  Two  such  values  are  required  in  the  cal- 
culation of  the  correction  in  Subroutine  CNVG.  Program 
control  then  returns  to  PS3,  where  now  ICVG  = 0 and  IPR  =0, 
and  variables  are  initialized  at  their  launch  values.  The 
program  is  now  ready  to  solve  the  second  problem. 


-21- 


TABLE  II 


Data  Required  for  Problem  1 in  Text 


DATA  SUPPLIED 

STATEMENT  NO 

. (#) 

COMMENTS 

CA,  TLZ , TLD 

3500,  3505, 

3510 

CA<0 

42012 

3626 

JCV=IMA-JNS=0 

CA, TLZ, TLD 

3500,  3505, 

3510 

A new  CA 

42022 

3626 

0. 

3882 

Final  ^ 

0.000002 

3885 

Sample  toler- 
ance 

TABLE  III 

Data  Required  for  Problem  2 in  the  Text 


DATA  SUPPLIED 

STATEMENT  NO 

'.  (#) 

COMMENTS 

CA,  TLZ, 

TLD 

3500,  3505, 

3510 

First  stage 
values 

24002 

3626 

CA,  TLZ, 

TLD 

3500,  3505, 

3510 

Third  stage 
values 

24212 

3626 

CA,  TLZ, 

TLD 

3500,  3505, 

3510 

A new  CA  - 
third  stage 

24222 

3626 

0. 

3882 

0.000002 

3885 

22 


Problem  2 : 


With  a specified  tilt  in  the  first  stage, 
the  tilt  in  the  third  stage  must  be  adjusted  so  that  the 
final  orbit  injection  heading  is  zero.  Table  3 specifies 
the  data  for  this  problem.  The  first  stage  tilt  value  is 
read-in,  the  first  two  stages  are  integrated,  with  the 
trajectory  values  being  printed  out,  and  the  final  re- 
sults for  the  variables  are  stored  for  use  as  initial 
values  in  the  integration  of  subsequent  stages.  'hen 
trial  tilt  parameters  for  the  third  stage  are  read-in, 
and  the  injection  conditions  are  calculated  and  printed 
out  as  a result  of  the  integration  of  the  third  and  fourth 
stages.  Based  on  this,  a new  CA  for  the  third  stage  is 
entered  by  the  user  along  with  the  other  tilt  parameters. 
Then  the  program  is  put  into  its  se If -correct ing  mode, 
iteration  on  CA  for  the  third  stage  ensues,  and  finally 
the  converged  launch  trajectory  is  printed  out.  Control 
is  then  returned  to  PS3,  where  the  variables  are  initia- 
lized at  their  launch  values.  The  program  is  now  ready  to 
solve  the  third  problem. 

Problem  3 ; With  specified  tilt  functions  in  the  first  and 
third  rocket  stages  and  a specified  coast  time  for  the 
second  stage  (TBD(2)),  it  is  desired  to  adjust  the  coast 
time  of  the  fourth  stage  to  obtain  ^ = 0 as  the  orbit  in- 
jection heading.  The  data  supplied  by  the  user  is  shown 
in  Table  4.  By  putting  JCV  = 2,  the  correction  procedure 
is  now  on  TBD (4) . It  proceeds  analagously  to  that  on  CA. 
After  converging  to  the  correct  coast  time  for  the  fourth 
stage,  program  control  reverts  back  to  PS3,  where  variables 
are  initialized  at  their  launch  values  again,  and  the 
fourth  problem  is  now  ready  to  be  solved. 

Problem  4 : With  specified  tilt  functions  in  the  first 
and  third  rocket  stages,  it  is  desired  to  adjust  the  coast 
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TABLE  IV 


Data  Required  for  Problem  3 in  Text 


DATA  SUPPLIED 


STATEMENT  NO. 


COMMENTS 


CA,  TLZ , TLD 

3500,  3505,  3510 

First  stage 

24012 

3626 

CA,  TLZ,  TLD 

3500,  3505,  3510 

Third  stage 

14212 

3626 

0. 

3500 

These  tilt 

0. 

3505 

parameters 

play 

0. 

3510 

no  role 
(no  thrust) 

143122 

3626 

TLS 

8687 

New  coast  time 
fourth  stage 

0.,  0.,  0. 

3500,  3505,  3510 

143222 

3626 

0. 

3882 

.000002 

3885 

TABLE  V 


Data  Supplied  by  User  for  Problem  4 in  Text 


DATA  SUPPLIED 

STATEMENT  NO.  (#) 

COMMENTS 

CA,  TLZ , TLD 

3500,  3505,  3510 

First  stage 

14001 

3626 

o 

• 

% 

o 

• 

% 

o 

• 

3500,  3505,  3510 

Second  stage 

14114201 

3626 

Coast  time 
iteration 

2 

3845 

Temporary  NS 

TLS 

8687 

New  coast  time 

0.  ,0. ,0. 

3500,  3505,  3510 

Irrelevant  tilt 
again 

141242 

3626 

Sel^-correct 
coast  time 

PSSS 

3882 

4>  2 

.000002 

3885 

Tolerance 

Ca , TLZ , TLD 

3500,  3505,  3510 

Third  stage 

14204001 

3626 

4 

3845 

Restore  ori- 

ginal NS 

o 

• 

% 

o 

• 

% 

o 

• 

3500,  3505,  3510 

Fourth  stage 

143122 

3626 

TLS 

8687 

New  coast  time 

o 

• 

% 

o 

• 

% 

o 

• 

3500,  3505,  3510 

Irrelevant  tilt 
again 

143222 

3626 

0. 

3882 

Orbit  injec- 
tion heading 

.000002 

3885 

time  of  the  second  stage  to  obtain  a particular  heading 
il>2  at  the  end  of  the  second  stage,  and  then  to  adjust  the 
coast  time  of  the  fourth  stage  to  obtain  an  orbit  in- 
jection heading  ^=0 . The  data  supplied  by  the  user  is 
shown  in  Table  5,  The  rocket  launch  trajectory  is  printed 
out  for  all  four  stages  in  the  course  of  the  program. 

As  demonstrated  by  the  examples  discussed,  the  pro- 
gram is  flexible  enough  to  handle  a wide  variety  of  rocket 

problems.  Whatever  rocket  problem  the  user  is  dealing 

. 1 . 

with,  it  seems  advisable  to  check  the  data  to  be  supplied 
to  the  program  by  tracing  its  path  through  the  flow  chart. 


Example;  A Hypothetical  2-stage  Launch  Vehicle 

For  typical  input  data  the  numerical  output  from  a 
specific  example  is  obtained  next,  and  its  use  in  computa- 
tions is  demonstrated.  As  in  the  preceding  section,  the 
rocket  in  this  example  has  four  stages  of  motion,  the  first 
and  third  being  powered  and  the  second  and  fourth  being 
coasting  stages.  One  might  typically  be  given  the  rocket 
data  of  Table  6.  It  is  a simple  matter  to  calculate  the 
input  parameters  for  the  program  from  Table  6 with  the 
help  of  the  following  equations : 


SMZ 

SMD 

TH  ; 


T 

vac 


/I 


sp 


T 

vac 


VSLR  = x 
(TBD)  (SMD)  = Mp 


(24) 


The  left-hand  sides  of  these  equations  are  written  in  the 
notation  of  the  computer  program  (Appendix  B) , while  the 
right-hand  sides  are  written  in  the  notation  of  Table  6. 
In  Appendix  D the  solutions  of  Problem  1 and  Problem  2 
of  the  preceding  section  are  obtained  on  a remote  time- 
sharing DEC-10  terminal  for  the  hypothetical  rocket  of 
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Rocket  Data  for  a 

PARAMETER 

TABLE  VI 

Hypothetical 

STAGE 

2-Stage  Launch 

FIRST 

Vehicle 

SECOND 

Ignition  Wt.*MQ 

(ibs) 

400000 

120000 

Stage  Wt.  *M 

(lbs) 

274000 

110000 

Total  Propellant 

Wt.  a M 

D 

(lbs) 

264000 

105000 

* 

Vac.  Spec.  Impulse 
* Jsp 

(sec) 

302.3 

314.3 

Vac,  Thrust  s T 

vac 

(lbs) 

665000 

220000 

Vac. Sea  level 

Thrust  Ratio  s x 

Stage  diameter 

(ft) 

1.15 

10 

10 
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* 


Table  6.  The  user  enters  data  for  the  rocket  problem  when 

the  program  prompts  him  to  do  so  with  an  ENTER message. 

The  value  10  seconds  is  arbitrarily  used  for  the  coast  times 
of  the  second  and  fourth  stages  of  the  rocket  motion.  This 
might,  for  example,  allow  enough  time  between  burnout  of 
one  powered  stage  and  ignition  of  the  next  powered  stage  to 
take  care  of  jettisoning  launch  vehicle  vestiges,  in  this 
particular  exercisethe  user  is  trying  to  obtain  the  orbit  in- 
jection velocity  at  zero  heading  (in  this  case,  the  perigee 
of  the  orbit)  and  at  an  altitude  of  100  nmi.  = 185,2  Am. 

In  order  to  find  this, the  user  plots  the  three  output  values 
altitude  vs.  velocity,  as  in  Fig.  7 (x  - marks),  and  then 
interpolates  (to  o— mark)  to  find  an  injection  velocity 
v*  m 7.833  Km/sec.  Incidentally,  the  total  computer  runtime 
for  the  output  in  Appendix  D was  26  seconds. 

The  value  just  found  for  the  orbit  injection  velo- 
city at  100  nmi.  altitude  is  in  a stationary  earth  approxi- 
mation, as  discussed  in  the  Introduction.  The  earth's  rota- 
tion is  easy  to  account  for  C 1]  in  separate  computations . To 
do  this,  one  notes  first  that  the  total  ground  range  covered 
in  this  example  is  about  670km. , which  is  small  enough  that  a 
flat  earth  approximation[l]  is  in  order.  Then  the  actual 
orbit  injection  conditions  for  the  rocket,  including  earth ' s 
rotation,  are  given  by 

v - [v  ^ + 2v'v  sin  va  + v 
o o o 


0 - 0° 

tan  a = vsin  ao  + vo 


v cos  a 

o 

where 

vo  s wer  Cl  + fy/R]  cos  L 


(25) 


Here,  v is  the  injection  velocity  in  the  stationary  earth 
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approximation,  vQ  is  the  effective  initial  velocity  in  an 
easterly  direction  imparted  by  earth's  rotation,  aQ  is  the 
launch  azimuth  angle  as  seen  by  an  earth  observer , which  is 
established  shortly  after  launch  by  suitable  pitch  and  yaw 
rotations,  and  a is  the  actual  azimuthal  angle  at  orbit 
injection  as  seen  by  an  inertial  observer.  These  azimuthal 
angles  are  measured  clockwise  from  the  local  north  direction 
to  the  direction  of  motion  projected  on  the  local  hori- 
zontal plane.  In  the  expression  for  v , u>E  is  the  angular 
rotation  rate  of  the  earth,  R is  its  radius,  y^  is  the 
rocket  altitude  after  first  stage  burnout,  and  Lq  is  the 
latitude  of  the  launch  site.  The  orbit  inclination  angle 
i,  which  is  the  angle  made  by  the  normal  to  the  orbital 
plane  (determined  from  the  right-hand  rule  by  curling  one's 
fingers  in  the  direction  of  motion  of  the  satellite)  with 
the  earth's  rotation  axis,  is  determined  from  the  pre- 
ceding parameters  by 

cos  i = sin  a cos  (26) 

O 9 

where 

o * a s 2tt  and  Jlq]  s is  if  - | Lq  | 
in  radian  measure.  By  way  of  illustration  it  is  supposed  that 
the  rocket  in  the  example  is  launched  from  Cape  Kennedy, 

for  which  L = 28.5°.  One  has  that  v = 7.833Km/sec,  cu_  * 
o & 

7.292116  rad/sec.,  R = 6378. 145km., y1  « 64/:m. , and  aQ  is 
variable  in  the  example.  Fig.  8 is  a plot  of  a and  i vs.  aQ 
and  in  Fig.  9 is  a plot  of  v vs.  i.  This  is  as  much  informa- 
tion as  is  needed  for  the  launch  problem  example.  It  will  be 
seen  from  these  figures  and  Fig.  3 that,  as  expected,  the 
apogee  altitude  will  be  maximum  when  the  launch  direction 
is  due  east  (aQ  = 90°) ; then  it  is  about  1075nmi.  On  the 
other  hand,  for  launch  directions  in  the  range  190°  s aQ  s 
350°  there  is  a question  of  the  stability  of  the  parking 
orbit,  since  in  this  range  orbit  altitudes  dip  below  lOOnmi. 
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DISCUSSION 


The  program  developed  here  was  employed  in  calculations 
on  large  launch  vehicles.  Curves  of  the  type  of  Fig.  7 
were  found  for  the  injection  conditions  in  the  stationary 
earth  approximation.  Earth's  rotation  effects  were  sub- 
sequently included  in  the  manner  described  previously. 

The  results  were  compared  to  corresponding  calculations 
which  leave  out  explicit  consideration  of  thrust  increase 
and  drag  effects,  but  rather  assume  that  they  cancel  in  an 
approximation  delineated  by  Ruppe  [4].  A similar  approx- 
imation was  used  by  Ehricke  [ 3] . The  curves  based  on  this 
approximation  yielded  injection  speeds  at  a given  altitude 
2-3  % less  than  those  calculated  with  the  present  program. 
For  very  large  rockets  the  thrust  increase  effect  outweighs 
the  drag  effect  by  about  this  amount.  Because  of  the 
sensitive  dependence  of  apogee  altitude  on  perigee  injection 
speed  for  eccentric  orbits,  it  was  concluded  that  explicit 
inclusion  of  thrust  increase  and  drag  effects  was  necessary. 
It  was  also  possible  to  calculate  results  with  the  present 
program  for  a large  launch  vehicle  which  had  been  indepen- 
dently calculated  elsewhere.  These  were  results  of  the 
nature  of  Fig.  9,  and  the  agreement  between  the  two  sets 
of  calculations  was  within  0.05%.  While  this  kind  of 
agreement  is  probably  fortuitous,  it  is  also  somewhat 
encouraging. 

Realistically,  one  would  expect  the  present  program 
to  yield  information  about  large,  unfinned  rocket  systems, 
consistent  with  the  data  base  to  within  1%  for  the 
injection  speed.  The  neglect  of  angle  of  thrust  (j3  in 
Fig.  2)  is  supposed  to  be  a good  approximation  for  large 
rockets  [3,  4],  and  a particular  form  for  the  angle  of 
attack  profile  is  apparently  unimportant  for  these  systems. 

A variation  of  parameters  in  the  profile  used  here  produced 
~ 0.2%  errors.  Beyond  this,  one  can  argue  that  small  errors 


-30- 


in  the  aerodynamic  coefficients  or  atmospheric  model  are 
not  important,  because  the  total  contribution  of  aero- 
dynamic effects  is  relatively  small  for  large  rockets  [ 4] . 
Another  consideration  is  that  significant  improvements  on 
the  present  program  are  probably  not  possible  for  many 
rocket  problems.  Data  for  such  a problem  are  often 
comparable  with  that  of  Table  6.  Further  refinements  on 
the  approximations  would  require  more  data  than  is  given. 

One  would  need,  for  example,  drag  and  lift  coefficients 
specific  to  the  rocket  system,  details  of  the  angle  of 
attack  or  thrust  angle  profiles,  the  time  history  of  thrust 
and  mass  variations,  and  data  on  the  atmospheric  conditions. 
Furthermore,  uncertainties  in  rocket  parameters  in  the 
data  base  can  and  often  do  amount  to  a few  percent,  and 
this  is  therefore  an  inherent  limitation  to  accuracy. 

The  present  program  is  certainly  flexible  enough  to 
handle  a large  variety  of  rocket  problems.  From  the 
preceding  discussion  it  appears  to  be  entirely  satisfactory 
for  most  launch  system  analyses,  particularly  those  which 
involve  large,  unfinned  rockets. 
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Fig.  2 - Force  diagram  for  the  rocket  with  symbol  legend 
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READ  TH(II,  I - 1.NS  (VAC  THRUST- LBSI 

READ  VSLRID.  I - l.NS  I VAC  TO  SEA  LEVEL  THRUST  RATIOI 

READ  TBDIII.  I - l.NS  (BURN  DURATION  TIME  FOR  EACH  STAGE-SECI 

READ  Dl(l),  I - l.NS  (DIAMETER  OF  STAGE-METERSI 


69  READ  Nld),  I - l.NS  (NO.  OF  INTEGRATION  INTERVALS-TIME  MESH) 


COMPUTE  TBSd),  I - l.NS  (TIME  AT  BEGINNING  OF  STAGE-SEC) 


3 INITIALIZE  PARAMETERS  AS  APPROPRIATE  TO  LAUNCH  CONDITIONS 


'CVG  - 1/TIs 


ICVG  « 2 


IS 

ICVG>  1 

? 


REAO  ANGLE  OF  ATTACK  (a)  PARAMETERS  CA  IAMPUTUDE-RADIANS), 
TLZ  (TIME  AT  START  OF  a#0-SEC>,  AND  TLD  (DURATION  OF  a#0-SEC) 


READ  PARAMETERS  WHICH  SPECIFY  EXECUTION  OPTIONS:  NT,  I0P,  ICE. 

IBO  IBP  IP\  / III  A I AlO 


COMMENCE  NUMERICAL  INTEGRATION  IRUNGE-KUTTA  METHOD)  OF  THE  EQUATIONS 
OF  MOTION  FOR  NT  ROCKET  STAGES 


Fig.  6(a)  - Flow  chart  for  RALLY  6 (Dec- 10  version) 
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SOME  NUMERICAL  integration  details 
DO  10  II  - I NT  [SETS  UP  DO  LOOP  TO  INTEGRATE  NT  ROCKET  STAGES) 
PARAMETERS  ARE  SET  UP  FOR  STAGE  I I - ICE  * III.  ANO  ALL  ENTITIES  ARE 
CONVERTED  TO  MKS  VALUES.  SET  UK  - I 

00  20  J - 1.NN  INN  TIME  INTERVALS  OF  WIOTH  H SEC  ARE  INTEGRATED  OVER  IN 
THE  SOLUTION  OF  THE  EO  NS  OF  MOTIONI 
PRINT  OUT  VARIABLES  AT  EACH  J IF  IPR  - 0 
20  CONTINUE 


10  CONTINUE 


CALL  STORE  (WHICH  STORES  RESULTS  OF  THE  LAST  INTEGRATION  STEP  IF 
ICD  - 1 OR  INITIALIZES  THE  PARAMETERS  AGAIN  TO  THEIR  VALUES 
PRECEDING  THE  INTEGRATION  OF  NT  STAGES  IF  ICO  * 21 


SET  YN  - PSF-PSSS 


-YES 

XN  - CA 


XN  - TBO  INS) 


CALL  CNVG  (WHICH  COMPUTES  A CORRECTION  TO  XN  BASED  ON  XO  XN 
YO  AND  YN  THEN  SETS  XO  - XN.  YO  - YN  IF  THE  CORRECTION  IS  LESS  IN 
MAGNITUDE  THAN  EPS.  THEN  IPR  - 0 AND  ICVG  * I ARE  SET 
A NEW  XN  IS  COMPUTED  I 


Fig.  6(b)  - Flow  chart  for  RALLY  6 (Dec-10  version) 
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a (deg) 


a0  (deg) 


Fig.  8 - Injection  azimuth  angle  and  orbit  inclination  angle  dependence 
on  launch  azimuthal  angle  for  the  example  rocket  (see  text) 


i(deg) 


v (km/sec) 


APPENDIX  A 


Analytic  Atmospheric  Density  and  Pressure  Profiles 

In  hydrostatic  equilibrium  we  have 

dp/dy=  - pg  , (Al) 

where  p is  the  pressure,  y is  the  altitude,  p is  the 
density,  and  g is  the  acceleration  due  to  gravity.  The 
ideal  gas  law,  Eq.  (13)  of  the  text,  is  reproduced  here  as 

p = pRT/M  , (A2) 

from  which  Eq.  (Al)  is  rewritten  as 

dp/p  = -dy/H  (A3) 


where 

H = RT/Mg 

is  the  "scale  height".  From  Eq's  (A2)  and 
dp/p  = dp/p  - dT/T  + dM/M  = dp/p  - dH/H 


(A4) 

(A3) , we  get 
- dg/g  (A5) 


There  are  several  numerical  tabulations  of  standard  at- 
mosphere density  and  pressure  profiles,  and  if  we  choose 
one  of  them  [6],  then  for  an  altitude  y near  a particular 
tabular  value  y , we  have  from  Eq.  (A4) 


H 


H 


o 


+ 0 

o 


(y-yQ) 


(A6) 


where  Hq  and  0 are  constants  evaluated  at  y , and  Eq. 
(A6)  is°simply°a  first  order  expansion  of  H §bout  y . 
Direct  evaluation  yields  ° 


1 

A 

_ " 

dH 

RT 

. dT  1 dM  . 2 

1 - + 

dv 

0 

>i 

ji 



Mg 

T dZ  M dZ  R+y 

(A7) 

o 


where  R is  the  radius  of  the  earth,  and  we  have  used  the 
known  variation  of  g with  altitude.  The  case  of  constant 
gradient  of  scale  height  is  a well  known  simplification 
in  the  integration  of  Eq‘^  (A4)  and  (A5) , which  become 


dm  _ -1  dH  p 

e “?0  ~h  “9  i~ 


(AS) 


■44- 


aua).  -/Vi\  “ fig_  -/JM 

P ^ pq  \ 0O  / H ' 'Vo  \ Ho  / 

(A9> 


These  equations  are  thus  interpolation  formulae,  which 
when  used  in  conjunction  with  some  of  the  values  in  the 
numerical  tabulation  [6],  become  the  desired  analytic  ex- 
pression of  atmospheric  density  and  pressure  (see  text) . 
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APPENDIX  B 

DEC-10  Version  of  RALLY  6 


. TYPE 
001  00 
00110 
00120 
00300 
/ 

00400 

00900 

07. 

00*00 

00700 

00300 

5-04, 

00300 

01000 

239* 

01  too 
01200 
35-03. 
01300 
/ 

01390 

01379 

01400 

01900 

01990 

01*00 

01*29 

01*37 

01*90 

01*79 

01700 

01729 

01790 

01300 

01390 

01379 

01300 

01390 

01379 

02000 

02009 
0200* 

02010 
02090 

02  079 
02100 
02190 
02179 
02200 
02290 
02279 
02300 
02  379 
02400 
02900 
02*00 
02*09 
02*  1 0 
02*19 
02*20 
02*29 


RALLY*. FOR 

COMMON  '3 1 '*5C * 3.  S3*  VSR*  RHO . PRZ . P . VT 

COMMON- Bl'YS  <3*  * YD  *3  • .9HZ<3*  .ALZ<3>  .MZ‘3<  .BT<3*  .RR<3* 
COMMON 'Cl  SMZ<*»  * SMD  >*» * TBit*- * TBD.** • Dt  <** • TH  <*  • « NI 
DATA  < YD  <I»*l*l*3» '10. *29. .49. .99. *39. *30. .100. . 120. . 140. 

DATA  < YZ  * I ■ . 1*1. 3* '10. . 10. .29. -99. * 39. * 39. *30. . 100. . 120. 
DATA<RMZ<I • . 1*1.3* -'4. 17*5-04.4. 17*5-04.4. 0435-09. 9. *505- 

13. 1335-03.3. 1335-03*3. 3425-03. *. *425- 1 0. 3. *1 35- 1 1 ' 

DATA  <ALZ < I .• . 1*1.3.  '-3.  13035-04.-3.  13035-04.-3.  12395-04. 
1-3. 10335-04.-3. 03495-04.-3. 03495-04. -3. 03215-04.-3. 0373 


120 


130 

140 

192 


190 


1*0 

*3 

170 

130 
1 00 
200 
300 


2-3. 07735-04' 

DATA <HZ < 1 • * 1*1.3.  '*.  4037.*.  4037. *.4333.  7,3707.9.7239.9.7 
19.7323**. 3437*10. 733*' 

DATA* FT  < I . , 1*1.3* '-2. 03735-01.2. 00*45-03* *. 3345-02. 2.  44* 

1-7.22975-02* 1.77115-03.*. 144*5-02. 2. *32*5-01 . 2. 33905-01 

DATS  <99 < I • . 1*1 *3*'2.*1 4509. 2. *145 09. 2.9345 04. 4. 2335 02. 
19. 02.9. 02*2. 033*4. 0095-01* 4. 3245-02' 

SC* 3. 33*012503 
9**373. 149 
Pr*3. 1 41932*34 
°RZ*1 . 01 3233505 
ICV'3»0 
IP9*0 
TYPE  120 

FORMAT  <’  ENTER  NS'. 

READ  <3*100.  NS 
MM S *NS 
TYPE  130. NS 

READ  '3*200.  < SMZ<I.»  * 1*1, NS» 

TYPE  140. NS 

cORMAT  < ' ENTER  * II  * VAL'JSS 
READ  <3*200.  < SND < I . * 1*1 . NS., 

TYPE  150. NS 

FORMAT  <'  ENTER  '*11.'  VAL'JSS 
READ  <3*200.*  <TH<  I * . 1*1  *NS> 

TYPE  132 

FORMAT  <'  ENTER  VSLR*. 

READ  <3*200*  VSLR 
TYPE  1*0* NS 

FORMAT  ('  ENTER  '.II.  VAL'JSS 
READ  <3*200.  < TBD< I • • 1*1 .MS. 

TYPE  170. NS 

FORMAT  ('  ENTER  .11*  VAL’JSS 
READ  <5*200.  <DI <1  - . 1*1 .MS. 

TYPE  130* NS 

FORMAT  <'  ENTER  '*11.'  VAL'JSS 
READ  <3.300*  <NI <I> * 1*1. NS. 

FORMAT  <'  ENTER  '*11.'  VAL’JSS 
cORMAT  <11* 

FORMAT  <F 1 0. 0 • 

FORMAT  <15* 

;*0. 

T B * < l .*  * 0 . 

Ie  <NS.SO. 1*  30  TO  3 
DO  3 1*2. NS 
<*X+TBD< 1-1 » 


FOR  SMZ' • 

cOR  SND' ■ 

FOR  VAC.  TH' . 

FOR  TBD'  • 

cOR  D I • • 

FOR  NI'. 
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03s  30 

TBS  < I » -X 

03700 

3 

YI-0. 

037?  0 

TI-0. 

0377? 

VI-O. 

03737 

XI -0. 

03300 

PSI-9I  -3. 

03300 

IP-l 

03300 

y < V SLR- 1 .<  - V SLR 

03*3? 

9L-0. 

03*?0 

3? 

19  (ICV'S.SQ.  1;  ICV-3  - 3 

034?3 

Ic  <ICVS.  ‘3E.  1*  -53  TO  31 

034?? 

19  < 1 99 . -3T . 1 < -50  TO  310 

03430 

TYPE  190 

0347? 

190 

90RH9T  <'  EHTER  VSL’JES  ! 

03?  00 

RE9P  <?«  300'  ‘29 

0 3?  0? 

9E9D  <?«  300'  TLZ 

03?  10 

9E9P  <?«  300»  TLB 

03??  0 

TY9E  33? 

03?7? 

33? 

90RH9T  <'  EH  TER  VHL'JES  ! 

03?73 

1 J!1  3 ' < 

03333 

RE9D  <?.?00»  HT.IOP.ICE 

03300 

?00 

90RH9T  <311  • 

03-31? 

19  < JUS. £9. O'  50  TO  303 

03330 

TY9E  130 

0334? 

RE9B  <?,100»  MS 

03330 

700 

90RM9TC  I NT.  T I ME < SEC 

03370 

l.>  -39.  RM-3 <K!1>  ' ' 

03373 

303 

19  < 199. LT . 3.'  50  TO  310 

03373 

TYRE  13? 

03373 

13? 

90PM9T  <'  EHTER  VHl'JSS  ! 

03333 

9S9D  <?*  300 • PS  S3 

0333? 

RE9B  <?. 300 • E9S 

il  33-43 

YO-PSP-®CSS 

03334 

310 

COMTIM'JS 

03300 

DO  10  II-l.MT 

04000 

I-ICS-I I 

040?0 

IJK*  I 

04100 

DP I -HI < I • 

04300 

H-TBBCD-'OBI 

043?  0 

VT-TH  < I 1 -4.4433 

04330 

S-PI-BI  (l.i  -PI  <1 ' <4. 

04300 

HH-HI  <I> 

04  3 0? 

9'1-SH2< I " 3.30433 

043?3 

1 1 1 - 0 

04  333 

19  <199.  ST.  1;  -50  TO  413 

04  337 

TYRE  700 

04373 

413 

COMT  IH'JE 

04400 

DO  30  J-l.HH 

044  03 

III-J 

04404 

19  <VI.»  33.31.33 

04  403 

31 

T-VT'VSLR 

04403 

9P-0. 

044 1 0 

9 J 1 — 0 . 

04413 

YL1-0. 

04414 

XH 1*0. 

04*13 

C SLL  V I NCR  < P S I • T • R V I • 1 

0441-3 

VK 1 -H-VK 

04430 

-30  TO  3 3 

04433 

55 

■:*ll  DEH < ID* YI; 

04434 

C9LL  c0RCES<VI.9l»S*T.9 

04433 

C9LL  VIHCR  <PS I * T*  9:1.  Y I . 

04440 

C9LL  PS IMCR  <VI . 9S I • T. 9M 

9LT.  '1 


51 


04  4 4? 

:pll  plt  <vi.psi.yy> 

04447 

CPLL  'jPP'j'VI.PSI.Yl.X) 

04443 

VK 1 »M*VX 

04450 

P J1 *m*?sS 

04  455 

YL 1 »M4YY 

04477 

XM  1 »M*X 

04500 

23 

TB»TI*Hx2. 

04-300 

V l »v  I ♦vs  12. 

04700 

PS1*°SI*PJ1  2. 

04  300 

Y 1 »Y I *YL 1^2. 

04  300 

■!PLL  D£M<  ID. Y l > 

05  0 0 0 

■:pll  tilt  <•:*•  tl::.  tld.  tb.pl  ■ 

05100 

SPLL  MPSS  <T?.  I • PM* 

051  03 

SPLL  F39CSS<V1. PL. S.T. PDSL- 

05200 

'S  P'-L  VIMCPYPSl.T.pn.Yl.FD.ftL. 

05300 

CPLL  P 1 1 MS  9 < V l • PS  1 . T . PM.  t 1 . F i. 

05  4 0 0 

S ML  L PLT  « VI  » °S  1 * YY  * 

05450 

IPLL  599'5  ■ V 1 • P 3 1 • Y 1 . X • 

05500 

V‘:2*H*V,> 

05-300 

° J2*H*CSS 

05700 

yl2*m*vy 

05750 

-!M2*M*X 

05300 

V2«VI ♦VP 2- 2. 

05300 

°'2*P? I*P J2  '2. 

0*3  0 0 0 

i'2*i'I*YL2  2. 

03  1 0 0 

CPuL  DEM  < I D«  Y2» 

0-3 1 02 

CPLL  F39CES  <V2»PL»  '. T.FD.FL  - 

0-32  0 0 

CPLL  VIMCP  <P32*  T.  PM. Y2.FD.Pl. 

0-3300 

: PLL  PS  I MCP  1 V2 . P S2  . T . PM.  Y2  • FL 

03  4 0 0 

cpll  PLT  • V2. P S2«  YY • 

03450 

CPLL  -5PP'5  * V2«PC2»  Y2*  X • 

03500 

VY 

033  00 

= J3*M*°SS 

03  7 0 0 

YL3»H*YY 

03  75  0 

\'M3«H*X 

03300 

V3*VIW>-  3 

03  3 0 0 

PS3*P5I*P  '3 

07  000 

Y3«Y  I ♦YL  3 

071  00 

TB»T I 

07200 

CPLL  DEM  Y ID«  Y 3 • 

07250 

CPLL  TILT  YCP.TL2.TLD.TB.PL> 

07275 

CPLL  MP3S  <TB- I.  PM' 

07277 

1 P'_L  F3PCES  '*  V 3* ML . 3.  T.FD.FL.1 

07300 

CPLL  V I MC 9 * P S 3 * T.PM. Y3.FD.PL. 

07  400 

S PLL  c' S I MCP  < V 3 . P$  3. T.  PM.  Y 3 . FL 

07500 

SPLL  PLT  YV3.PS3.YY* 

07502 

CPLL  -5PPS  < V3>  PS 3.  Y3»  X.» 

07505 

VK4*H«VK 

07510 

?J4*M*PS3 

07515 

YL  4*M4YY 

07557 

XM4*M*X 

073  00 

VF *V I ♦ 1 VK  1 >2.  ♦ ' VS  2 ♦VS  3 1 ♦VS  4 * - 

07700 

? ! F *°S  I ♦ < P J 1 ♦ 2 . • • P J2  »P  J 3 • ♦?  J 4 

07300 

YF «Y I ♦ • YL 1 ^2 . * < YL2+YL 3 • ♦ YL  4 • <* 

07325 

IF  «X I ♦ • XM 1 +2 . • < XM2+XM 3 1 *XM4 • 

07350 

IF  ' I °P.  '5T . 01  33  T3  501 

07300 

TYPE  300. J. TB.VF.PSF.YF.XF 

07350 

*3»00 

C3PMP  T y 1 X . 1 5 . 5 ' 1 X . E 1 0 . 4 > < 

07355 

501 

TI*TB 

07330 

V I »VC 

07335 

?SI*°CF 

07370 

YI-YF 

s » 


• HL • PS  s 


XHIS  PAQK  IS  BEST  QUALITY  SUAUt**** 
htaM  CuFY  fUNUSHB)  10  S®0  *— - 


-48- 


e 


r 

I; 


J 


0*32? 

03?  00 

03-300 

03700 

03300 

03300 

033? 0 1 

0337? 

1 0000  3 

10100 
10200  ? 
10300 
10400 
10402 
10?00 
10??  0 
10*00  ? 
10700 
10300 
10300  It 
1 03?  0 
11000 
11100  3 

11110 
11120 
11130 
11140 
1 1 ISO  4 
11200 
11300 
11302 
11400 
11402 
ll?00 
1 1 ? 02 
1 1*00 
11*02 
11*04 
1 1*0* 

11700  7 

11300 
11300 
12000 
12100 
12200 
12300 
12400 
13020 
13030 
13400 
1 3?  00 
1 3*00 
13700 
13300 
13300 
1 330? 

13310 

1331? 

13320  ? 

14000  7 

14100 
14200 


RETURN 

END 

SUBROUTINE  TILT  <C.3«  T’_Z*  T'_D»  TB»3L.- 
3!_-0. 

IF  <TB.5T.T!_Z.  50  TO  1 
30  TO  ? 

TT  -TLZ*T’_D 

IF  <TB.5T.TT-  50  TO  ? 

T-TB-TLZ 

3L-C3*T*  <1 . -T/TLD- 

RETURN 

END 

SUBROUTINE  DEN<ID»Y» 

COMMON/* l /5C » R » S3  > VSR . PMQ.  PRO . P,  VT 

COMMON'S 1/YZ <3 > »YD<3>  <RHZ(3  - .3LZ<3-  .HZ <3-  » BT  <3- »PR<3> 
IF  < ID. 5T. 3>  50  TO  10 
IF  <Y.LE. YD<ID- - 50  TO  3 
I D-  ID*  l 

IF  <ID.LE.3-  50  TO  ? 

RHO-O. 

R^O. 

50  TO  7 

IF  <ID.!_T.  2'  50  TO  4 
YY-YD<ID-1 > 

IF  <Y . 5T . YY > 50  TO  4 
ID- ID- l 
50  TO  3 
X-Y-YZ  <ID> 

T-BT  < ID ■ 

H-HZ- ID-+T4X 
C-HZ  < ID • 'H 

3-RHZ  < I D < / < 1 . +3LZ  < I D » •X.t 
RP-RR  < I D ■- 
B— ;♦♦<  <1 . ♦ T> 'T  > 

1 . ' t * 

RHO-R*B*l 000. 

°-3P*BP'l 0. 

S S-JORT  < 1 . 4*R/RMO> 

SS-SS' 1 000. 

RETURN 

END 

S UBPOUT I ’IE  RSI  NCR  ' V . PS  * T . 3M . Y . F!_  - RL  • P S l < 

COMMON/ 31  ' 5C - P • S S » V SR - RNO . PPZ « R • VT 

X-9-fY 

B-X*X 

®S S*5C Y <V+B> -V -X 
RSS--PSS*C02  -PS- 

CD-CL*T 

PS  S-PSS+CDrS  IN  <*!_.> / <3M*V«»1 000. • 

RETURN 

END 

SUBROUT I NE  V I NCR  <RS • T . *M . Y , FD . 3L • VD • 

COMMON  /*  1 <'5C  • R » SS  • V SR  • RNO  - RRZ  • P • V T 
X-R+Y 

X«X4X 

IF <3L. £0.  0.  ■ 50  TO  ? 

TT*T*C0S<3L.--FD 
50  TO  7 
TT-T-FD 

VD  * T T / < RM# 1 0 0 0 . > - - 5C ' X - ♦ S I N < P S 

RETURN 

END 


.„15  MM  IS  MSI 4MLIII  PMC11MW4 
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1 4 300 

S'JBPO'JT  IMS  m‘_T  <V»  P i * YD  • 

14400 

YD*V*SIN«PS* 

14500 

9£T'J9fi 

14*00 

EMD 

14700 

i'JBPO'JT  IMS  '5PP-5  < V*  PC.  y . X > 

14300 

COMMON '*1  -5C  « P • i i . ViP.PMO.PPO.P*  vt 

14300 

X*V*COS ‘PS  • 

l 4320 

xx*cos ‘Pi* 

15000 

X*X*P^  <p*Y> 

15100 

PST'JPM 

15200 

SMD 

15300 

S'JBPO'JT  IMS  CMV3  <XN«  XO*  YO.  YM»  EPS.  IPP.  ICV-3> 

15400 

DY«  CYM-Y3  * '<XN-XO' 

15500 

DX-— Y!1*-DY 

15*00 

XO*XM 

15700 

YO*YM 

15300 

XN*XN+DX 

15300 

:3*PBS  <XO-XN» 

1*000 

IF  ‘.3.  SS.SPS*  '30  TO  5 

1*100 

I?P*0 

1*101 

ICVS*1 

1*200 

5 

PEVJfHi 

1*300 

SMD 

1*310 

S'JBPO'JT  I MS  STOPS  <V  I.PSI.YI.XI.TI.VVI.PS  SI.  YY I 

* I ID. I CD 

1*315 

l * 

1*320 

30  TO  <3.5*  ICB 

1*330 

3 

VVI*VI 

1*340 

PSSI*PSI 

1*350 

YY I *Y I 

1*3*0 

XXI*XI 

1*3*5 

TTI*TI 

1*370 

I ID*ID 

1*3  30 

•30  TO  7 

1*3  30 

5- 

VI  .'/VI 

1*400 

PS  I *PSi I 

1*410 

Y I . YY I 

1*420 

X I -»xx  I 

1*425 

TI=TTI 

1*4  30 

I D=  1 1 D 

1*440 

1*450 

SMD 

1*550 

S'JBPO'JT IMS  FORCES <V«SL»  S»  T«FD»F'_.> 

1**50 

‘1 'HI  ••  *5‘I « 9»  Z-  S 9 V SP • PH 3 «•  °9Z • P • V T 

1*750 

T *VT*  • 1 . -VSR*PyPP2.» 

1*300 

IF  <RHO.ES.  0.  f '50  TO  5 

1*350 

XM*V'SS 

1*350 

3t*XM- 1 . 1 4 

17050 

IF  <"XM.'3T.  0.  3'3-3-3  3 * -30  TO  2 

17150 

B*-1*.*42-?*H*!3 

17250 

CD.. 25-*-. 377 3*£XP <B> 

17350 

•30  TO  3 

17450 

2 

IF  <XM.ST.  1.2'3*  '30  TO  4 

17550 

B*rlBi  ‘S' 

17*50 

B*B-**3 . 5 

17750 

CD*.  *273-30.  '302*B 

17350 

■30  TO  3 

17350 

4 

CD* . 1 73 1 3+ . 5 3 1 3 ' XM 

i 3050 

3 

IF  <X!1.'3T.  1.5*  '30  TO  * 

13150 

I L -4 , *5-  0 . 5*CQS  < 2 . 35r?^X!1  > 

13250 

I DL  - 1 • 52——  0 • 2'3'3*CD3  '2.3  32  3 1 ♦Xii  • 

13350 

30  TO  5 

13450 

•j 

CL*0,  3+*.  'XM 

51 


APPENDIX  C 

PDP-11/70  Version  of  RALLY  6 


COMBO*/ A I 'Of.P.SS.VSP.PMO.PBZ.P.VT 

CO«*n«/»i /y7(93,yo(B3,prz(B3,ai?(B),mz(B),pt(9),pr(9) 

COMM(lR/d /SMZ(P1 ,.ftMO(ft) .TPS (ft) ,TBn< A) ,01 (6) , TH ( 6 ) , N I ( 6 ) 

DATA  vn/1 o.  ,75. ,4ft. ,55. ,8ft. ,90., ton. , 120. , 140./ 

DATA  V7/10.  , >0.  ,?A.  ,ftft.  ,85.,Pft.  .90.  , 100.  , 120./ 

OAT*  BH7/4.  I 7AE-04,  A.  1 7AF-04 , 4 . 0 4BF-05  , 5 . 6 50E-07  , Q . I 9 3E-09 , 

IB.  1 O1F-0O,  3.B4?F-09,ft.64?F-10,  3.61  3C-1  1 / 

DATA  AI.Z/-7.1 308F-04.-3.1 3 08E-04 . - 3 . 1 2 35E-04 . - 3 . 1 08 9F-04 , - 3 . 0945C-04 , 
l-3.O94AF-O4,-3.O921F-O4,-3.0B73F-04,-3.0778F-04/ 

DATA  HZ /ft. 40P7, A. 4097,6.4388,7. 87.07,5. 7?  3ft, 5. 72  35, 5. 7323, 
16.34B7.IO.7BBA/ 

DATA  bt/-?.037ha-01 , 2. 00  ft 4 F-0 3, 6.9948-02,2. 44A9E-0 3.-7. 2257E-0 2. 

1 1 .771 1f-O3,A.144ftF-02,?.A82ftF-0l .2.335OF-01/ 

nATA  PP/2.A14FOft,2.ft1 4F05.?.534F04,4.283F02,5.02,5.02,2.098, 

1 4.OOSF-0J  ,4.  174F-0?/ 

Gr«3.BBA01 2FOB 
Psft 37B . 1 Aft 
PT*  3 . 1 4 1 ftO?hft4 
TTVsO 
IPBsO 

irvnso 

PB7.S1  .Ot  32ftftFOft 
TYPF  110 


110 

FORMAT  ( • OUTPUT  TO 
PFAD  (ft,|O03  TPP 

2 

TYPF  1?0 

120 

F0P“»T  ( ’ fmtfu  ns • l 
PFAO  (5,100)  vj 

100 

FORMAT  ( t 1 ) 

NN.ftaNS 

TYPF  1 3 0 , f S 

1 30 

FORMAT  ( > f«jTFR  • , T 1 
PFAD  (ft. 200)  ( S **  7 ( T ) 

200 

FOPMAT  ( f 1 0 , 0 ) 

TYPF  140, *S 

1 40 

Format  ( • f«it»  u i , I 1 
MF  AD  fft,?00)  (S-»p{  t ) 

TYPF  1ftn,*S 

ISO 

FOpy at  ( • F'.T*o  • , r 1 
PFAD  (ft, 2003  (TM(T>, 
TYPF  )ft2 

15? 

FORMAT  ( > FN  — F P V.ftt.P 
BFAO  (ft, 200)  Vftl P 
TYPF  ) ftd  , U,ft 

IftO 

FORMAT  (•  FMTFP  • , T 1 
PFAD  ( ft , 200  3 ( TP  0 ( T ) 

TYPF  170, MS 

170 

FOPMAT  (•  FMTFP  • , T 1 
PFAO  (5.2003  ( 0 T ( T 3 , 

69 

TYPF  1 U ,3  , *\S 

1 BO 

FORMAT  (•  f\tfi»  >,t1 
PFAO  (ft,  300  3 f M T ( T > . 

300 

FOP  M A T ( ] ft  3 

X *0  . 

TftS ( 1 )*0. 

IF  (1’S.FO.  13  TO  TO  3 

no  ft  t *? , ''ft 

X*X ♦Too ( (-1 3 

ft 

T**ft(T3s< 

3 

Y T *0 . 

TTsO. 

V T »0  . 

VAf.UFft  FOP  SMZ'1 


VALMFS  Frp  SMD'3 
:1 ,9S1 


v At, UPS  Fpo  v*r.  T«. ' ) 


, • V A t.llFS  F('R  TPD  1 3 
, T = 1 , MS) 


V A t,  1 1 F S FOB  DI  • ) 


VAt.UFS  FOP  *1 T 1 ) 


/ 


,y 


// 


-/ 
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> 


L j 


25 

1 90 


225 

500 


700 
8 0 3 
1 05 

210 


41  7 

21 


2? 


21 


XI«0. 

PSI«PT/7. 

Tt>»  1 

VSP»(  VSl.9-1  1 / VSl.P 
AL»0  . 

TF  ttovg.fo.i)  rcvr;«2 

IF  (TCVG.GF.l)  GO  TO  210 
IF  (TTY.GT.1  ) GO  TO  210 
TYPF  190 

FORMAT  ('  (TNTFW  VA|,IIFS  FOP  C A , Tt,2  . TLD  • ) 

PFAO  f 5 , 700 1 r» 

PFAO  15*2001  TLZ 
PFAO  ( 5 , 200  l TUI 
TYPF  225 

FORMAT  (•  FNTFP  VAI,1IC5  FOP  NT , IOP  , ICF , TTY  , IBG , JCV  , IMP  , JNS  ' ) 
PFAO  (5,5001  «T,TOP,TCF,ITY,IBfi.JCV,IP».JPS 

format  («rn 

TF  (J«S.F0.01  GO  TO  #03 
TYPF  170 
PPTNT  170 
PFAO  (5,100)  «S 

FORMAT  (7X.4HTMT..1X , 1 OHTT#C(SFC. ) , 1 X , 1 OMSPFFD ( K /S ) , 1 X , 

1 1 OHMF  AP  T *.'G  ( P ) ,?X,9HAr.T.  (FM.  ) , IX  , 1 OMGR . PNG ( F M ) ) 

TF  (TTY. IT. 2)  GO  TO  710 
TYPF  195 

FORMAT  (•  FMTFO  VAL'IFR  fop  PSSS.FPS’l 

PFAO  (5,7001  PSSS 

PFAD  (5,700)  FPS 

YOsPRF-PRSS 

COMTTN'IF 

DO  10  I T S 1 , NT 

I * TCF ♦ T T 

I.JFaT 

DOIaNT(T) 

HaTBn( T 1 /POT 
VT*YM( T 1*4 . A4B7 
S*PT*OT(T)»01(T)/4. 

NNaN  T ( T ) 

A"s5M7  ( T 1 /?.  704*.  7 
T T 1 = 0 

TF  ( TTY ,GT. 1 ) GO  TO  417 
TYPF  7ft0 

TF  (TPP.FO.il  PRINT  700 
GnNT  I N I'F 

no  70  ,t*i , w v 
TTI=J 

IF  (VT 1 77,71  ,77 

T=VT/V5LR 

F0  = 0. 

P.JlaO. 

Y1 1aO. 

XMlaO. 

CAM.  VT«.'GW(  RFT  ,T,  fM  , YT  ,FT),  AL,VK) 

VF 1 aM» VF 
GO  TO  77 

CALI  OrN  ( T P , Y T 1 

CAJ.T.  F>tuff  5 ( V T , A f.  , S , T , FT)  , F|,  ) 

CAM.  YTNGP(pRT,T,AM,VT  ,FP,A?,,VF) 

CAM,  05IPGR(VT  ,PRT  ,T,AM,YI  , FT, , Al. , PRS  ) 

CAM.  »IT  (VT.P5T.YY) 

C AM,  GROG  (VT.P5T.YT ,X1 
VF 1 =M»VF 
PJ1aH»PR5 
Yt 1 *M»YY 
X* 1 an*  y 
TP sTT  AM /? . 
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;H  IS 
• tfoM  OOi  ^ 


* " «v.  v onP  


p^dJUSHTO  TO  OOC 


VJsVT*VK J /?. 

PSI*PSI4PJ1 /2. 

Y13YT4YLI  /2. 

CALL  PFM(in.Yl) 

CALL  TII.T  fCA.TLZ.TLD.TB.AL) 

CALL  “ASSfTB.I.AM) 

CALL  FOPCFSCVl .AL.S.T.FO.FL) 

CALL  VThCPfPSI .T.AM.Y1 .FO.AL.VK) 

CALL  PSlNfRf VI ,PSt .T.AM.Yt .Ft .AL.PSS) 

CALL  ALT  ( VI  ,PSI  ,YY) 

CALL  CRPO  f V 1 ,PS1  , Y1 ,X) 

VK2=P»VH 

pj?sH*PS.«! 

Yt.2sM«YY 

X M?*4  * Y 
V?*VT*V*?/2. 

PS2aPST ♦PJ?/?. 

Y5s V T ♦ YI 2/2. 

CALL  PF«(Tn.Y2) 

cali.  rnpcpsf v?. al.s.t.fd.fd 

CALL  VT»iCO(Os7,t,jm,y?.FD,AL,VK  ) 

CALL  PST*'CPf  V2.PS2.T,  AM,  Y?,FL,  AL.PSS) 

CALL  ALT  f V2 , PS 2 , YY  ) 

CALI  CPPO  fV2,PS2,Y2,X) 

VK M * V K 

p,)3sH«PSS 

YL3=M*YY 

V3sVT*VY  3 

PS  3sP«T*ovIJ 

Y3sYT*YL1 

TPsTT+H 

CALL  nFwftO.YX) 

CALL  TTIT  (CA.FLZ.Tl.P.TB.AL) 

CALL  MASSfTB.T  ,A-n 
C ALL  Ff'PCF  S f V 3 . »L , S . T , Fn . FL ) 

C»Lt.  V TNCP ( PS  3 » T , AM , Y 3 » FO , AL  » VK ) 

CALL  PST'JCP(V3,PS3.Tf  am,  Y3.FL.At.  ,PSS) 

Ctrl.  AJT  tV3,pS3,YY) 

CAT.L  CRPO  ( V3.PS3.V3,  XI 

VF4sM«VX 

PJAsMAPK*! 

YL4=P»YY 

X*4=H*v 

VFsVY*  ('<*  \ ♦?.  * ( VK2*VK  3 )+VK4  ) / A. 

PSF*PS  T ♦ f o,M  ♦ 2 , * ( P.t ?*P.t  3 ) 4P.14)  /K  . 

YF*Yl4f YL142.*(YL?4YL3)4YL«)/F. 
XFsyT4(XMl«2.*rXM24YM3)4XM4)/A. 
tF  (TTY.OT.O)  fin  TO  SCI 
TYPF  AOO,,l,TH.vF,PSF,YF,XF 

IF  (TPP.fd.n  PPTt'T  B20,.1,TB,VF,PSF,YF,XF 
FORMAT  n X . IS, KM  X.Fin.4)  1 

TTsTP 
V T * V F 
psrsPSF 

YTsYF 

XI  = YF 

COMTIMLF 

TF  f Trip.  A F.  4 1 CO  TO  3 3 
IF  fTT.MF.»|T)  r.r  to  in 
I C I)  s 1 

TF  fT.f".*>S3  1 CPs  2 
TF  ( T CVC . fo.  2 ) TCOs 1 

CAM,  5 T o p r (VT ,PST .YT.XT.TT .VVT.PSST , 
1YVT.YYT.TT 1, 10,110, ICO) 

CO  TO  1 2 


I 


E § 

F4  X 


* S 
£ 


a 


T F ( T .C-  .*51  (tn  TO  VB 
CO"T|mif 

in«i 

T CzTCVC 
TCVOaO 

IF  fTC.OF.1)  OP  TO  31 
IF  fTTV.FO.O)  on  ‘"I  11 
TYPF  f.OO,TTl  ,T1.VF,PSr,YF#XF 

IF  (IPP.FO.11  PPTRT  POO.IIT.Tp.vF.PSF.YF.XF 

IF  (MY.FP.1)  OO  TO  31 

YNsPSF-PSSS 

IF  MCV.FO.I)  r,n  TO  VIS 

xn=c» 

00  TO  270 
x»sT«r>{i;s) 

CALL  C*'VGf  XN.XP,  Yn,Y*»,FPS.TTY.  ICVO) 

TYPF  JSp.xw.YM 

IF  fIPO.FO.1)  OPIKT  1SP.XK.YN 

FOPPAT  f XN  Ff).  '.Flt.H,1  , YO  FO.  ’.F10.4) 

TF  fjrv.co.V1  00  m VIS 

CAaXN 

00  TO  31 s 

TLSsXM 

On  Tn  a 1 7 

IF  fTf.F0.V1  on  TO  flis 

IF  fJfV.FO.V)  on  TO  110 

XOafA 

00  TO  317 

XOsTPOf  TJX ) 

TYPF  ISA 

FOPM  AT  f • CT>TF«  TLS  • ) 

PF  AO  fS.VOOl  TLS 
IF  (JCV.HF.V1  on  TO  11* 

Tf  (IJ«.Fn.H»jFl  on  rn  017 

TTTsTLS-TCOr  TJ* 1 

NAHSTJK*1 

nn  Rid  1s«AP,fJ*'S 

TflSf  T ) sTPS ( T ItTTT 

Tpnr  t.ia  1 sti  s 

if  fir. »t. 21  on  to  n7 

TOPsTRO 

IF  (TPA.FF.1  1 r,n  to  117 

TYPF  41S 

FOP  MAT  f*  F'tTFO  S«*TR'1 
RF  AO  f S , 700)  S*<  f« 
fin  417  Tsl  , VS 
SM2(TlsS»*7f TItSFTN 

On  TO  ( 30.3.S7, VS, VI  TOP 

CP«ITT*"'F 

fmo 

SIIPROUTIFF  MASS  (TP  , T , AM  ) 

fOMMn«/ri /SFZfM ,SMO(S) .THSfPI .TMOfft) ,DI (6) .THf M ,MI (6) 
SM*SY7(Tl-.«Mr)f  T 1*(  tP-TPS(I  1 ) 

AF«S»/V. 214*2 

RFTUPM 

ENn 

SfiPBO'JT  T*’F  TITT  ffl  ,TL2,Tr.n,TP,Al  ' 

ALSO. 

IF  fTP.CT.Tf  2)  rn  rn  i 
on  to  s 
tt*tlz*ti n 

IF  fTP.07.TT1  00  rn  s 
TsTP-TLZ 

Af.  3f  A * T • f 1 .-T/T|  01 

PET'IPM 

run 


. , • • diiaT  aiuui* 


siibooiitt'j*'  pf'  (tn.Yt 
rnMMOW/»  1 /GP,  W . SS,  VSB  ,PMf),  PB7  ,P,  VT 

CO«“ftN / « * /Y7.(  PT  . »noi  .BH7.CQ)  .HLZCP)  .HZf  9)  ,RT(9)  ,PR(B) 

TF  (TO.GT.P)  pn  TO  in 
IF  (Y.-tF.YPf  III  1 GO  TO  1 
TDaTO»  t 

f f (TO.»,F.9>  on  TO  S 
PPOsO. 

0 = 0. 

RtJ  TO  7 

TF  (TO.f.T.?)  GO  TO  4 
YYsYOf  Tr'-  1 ) 

Ip  (Y.GT.vv)  GO  TO  4 
T Pa  10“ I 
GO  Tn  7 
XaY-VZr  ?t>) 

TaPTf TP) 

HsMZf TO)*T»x 
C=UZ(  Tr’)/w 

4aRM7fI0  1/  n.*»l?(Tf))*X) 

APrPP  < T" 1 

AaC** ( ( 1 . ♦ T 1 /T ) 

PPaC* ♦ ( 1 . /T) 

PHPaA****!  OOP. 

Ps»P*PD/10. 

SSaSOPTf \ ,^*p/»un\ 

SSaSS/1 "00, 
rftiio” 

F*IO 

SUPPOUTTMF  PStNfP  fV,PS,T,»M,y,Ff  ,»|,,PS5T 
Cnxann/41 /GC.R.SS, VSR , RMO , PPZ , P , VT 
At'Gal  .57070M77 

XaP*Y 
Rs**» 

PSSaGC/( V*»)-V/X 
P$SX“PG.S*SIM(  4»'G“OS  ) 

rn  = F|,  + T 

pssspss*rr*si  oc  ) / r *»*v*  t ooo, ) 

RKTORN 
END 

5UMR01ITT  *i  F VTMGP  ( P S . T , AM , Y , FP , 4 1 ,VP) 
rO»»ON/» \ /or ,o , SS, VSO ,PMO ,POZ,P , VT 
4 "Ga  1 .S707o*V7 
X aB  ♦ Y 
Xa  X * X 

TF  (At... ‘0.0.1  GO  TO  * 

TTaT*ST'i(4‘lG-4l,)-FO 
GO  TO  7 
TTaT-FP 

VPaTT/f 4><*1 OOP. )“(GC/XT*STM(PS1 
BFTMPW 
FNp 

SHPPrmTT,',F  4 t.T  ( V , PS , YO ) 

YPa  V *S T V ( PS  ) 

OFT'IBn 

F»ip 

SII«POI'TT*'F  GPpG  ( V , OS  , Y , X ) 

roxxoii/ji  /or  ,p  , ss ,vso .rmo.pbz ,p , vt 

AWGa 1 .S7070A797 
X a V * S T 'i  ( 4'.G“PS1 
XaX*P/(P*Y) 

RFTMPM 
FMP 

SURBOlITT'if  C"VGf  XV  , XO,  YO,  YN.FPS,  TTY,  TfVG) 
PYa(YH“VOT/(XM-XOY  <C 

DXa-YV/PV  *0 


XHmXN 

vnxvN 

xnaxr'*nx 

AaAPS rtn-xfc ) 

if  r * .rs-.FDS)  on  rn  •> 

ITYaO 

icvcai 

RFTltPM 

r*n 

Sl’PPP'lTM  r STPCF (XT  ,PST,YT,XT,TT,VVI,PSSI, 
lYYT.XXI.TTI.TO.tTO.TCn) 

cp  tp  (j.s)  im 

VVIaVT 

PSStsPST 

YYTSYI 

XXIaYT 

TTTrTY 

XTPsTP 

OP  TP  7 

VTxVVI 

PSIaPSST 

YT*YYT 

XTxXXT 

TTaTTT 

TPaTTO 

PfTlIBN 

Fro 

SUBROUTINE  FPPPe-Sf  V.AU  ,5,T,FP,FL) 

Cnaanr  /» i /C,r , o , ps , vsp  , pho  , prz  , p , vt 
ANCal .570706377 
Ta  VT*  f 1 .-V«R«P/P9?) 

TK  (Rup.FO.n.)  on  TP  6 

XMxV/SS 

AsX“'-l . 1 X 

IK  fX«'.PT.o.PP«BOi  on  TO  2 
Pa-16. 6*?6*»*» 

CPa.?P+. 3773*FXP(P) 

CP  TP  3 

TK  ( X M . GT . i . 2P 1 CO  TO  4 
Pa APS  f * 3 
PsP* • 3.5 

CPa. 6777-30. 00?*P 
CD  TP  3 

CPa.  t 7<St  3*  .5  3t  P/X* 

IK  CX-.CT.1.S1  on  TO  6 

CL  = 4.6-0.S*STM  AM0-7.866*XP1 

CPI.al . «77-n.?OR*STN( »HO-2.BO?81*X») 

CO  TP  S 
Ct.*n.«**,./x« 

CP!.a-.?«  t P + 3.27 27 /»* 

CaBMP*V*v/7. 

C*C*S* 1 OOOOOO . 

FPaCP*CPr.*A|,*Ar, 

FP»FP*0 

FL*Ct,«r 

RFTIIQN 

FNP 


«*ST  QUALITY  PRAOIld*^ 
r>uT  c ppQK  IS  BBS  I tJUAii* 

• Bis  r*w  «tchio  10  DOC  — ■ 

COPY  iw  U* 


1 


APPENDIX  D 


.EXECUTE  PALLY*. FOR 
L I Ilk  l Ldsmn  j 

CLMKXCT  PALLY*  Execution! 
EHTER  AS 
4 

EATER  4 VALUES  FOP  SM2 
400000. 

13*000. 

120000. 

10000. 

EHTER  4 VALUES  FOP  SI1P 

2200. 

0. 

700. 

0. 

EHTER  4 VALUES  FOR  VAC.  TH 
**■5000. 

0. 

220000. 

0. 

EHTER  4 VALUES  FOP  VSLP 
1.  15 
1.  15 
1. 

1. 

EHTER  4 VALUES  FOR  TBP 

120. 

10. 

150. 

10. 

EHTER  4 VALUES  FOP  PI 
3. 


3. 


EHTER  4 VALUES  FOP  HI 
40 
10 
10 
10 

EHTER  VALUES  FOP  CA.TLZ.TLP 
-.  001 
10. 

25. 

EHTER  VALUES  FOP  HT* IQP. ICE. IPR. IBS. JCV. IMA.JHS 


42012 

1\ 

1MT. 

TIME 'SEC. > 

SPEEP  iK  'S!' 

HEAD IMS <R> 

ALT. *KM. > SP.PHS'KM' 

IHT. 

TIME  <SEC . > 

SPEED <KxS> 

HEADIHScR • 

ALT. <KM. > SR. RMS  < KM • 

a \ 

IHT. 

TIME  < SEC . > 

SPEED  UOS..' 

HEAPIHS ■ P • 

ALT.  CKI1.  • SP.  RHS  OCM 1 

5 \ 

IHT. 

TIME -SEC. * 

SPEED <KVS> 

HEAPIHS'R* 

alt.  <km. sp.phs<km> 

10 

. 2900E+03 

. *S22E+01 

. 12*IE«-01 

. S315E+03  . I*37E*03 

jj  ) 

EHTER 
-.  005 

VAr.UES  FOR 

CA.TLZ.TLP 

^ o 

10. 

25. 

EHTER 

VALUES  FOP 

HT, I OP. ICE 

. I PR. IBS*  JC 

V.  IMA.  J'HS 

3 £ 

42022 

ihj 

EHTER 

VALUES  FOP 

PSSS.EPS 

f~i  X1 

u. 

. 000002 

it 

10 

■ 2*?00E  + 0 

. 772AE+01 

. 3612E  + 00 

. .3373E+03  . *333E*03 

XH  EC. 

- . 0^  1 2 U*?E “ 02  9 YD  ! 

EC.  . 3*12£ 

♦ 0 0 

$jr 

3 

0.  ft 

CO  j. 

s a 

e e 
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I 3 


‘ i 


l 


•: 

i ' j 

1 


n 

■ 

4 

!.] 


10 

. 290 OS +03 

.31 17S+01 

. 1494E+00 

. 2 029E ♦ 03 

. 735 15^03 

£0. 

-. 773333 13E 

-02  . YD  EO.  . 1 494E+00 

10 

. 2900E+03 

.3330E+01 

. 31 00E-01 

. 1 143E+03 

. 7324E+0? 

£■3. 

-. 3033233 0E 

-02  . YD  £9.  . 3100E- 

01 

10 

. 2900E+03 

. 3417E+01 

. 3159E-02 

. 9207E+02 

. 7923E+03 

£0. 

-. 3 033323 OS 

-02  . YD  £0.  . 3159S- 

02 

10 

. 2900E+0  3 

. 3423E+0 1 

. 3940E-04 

. 395SE+02 

. 7940E+03 

£'?. 

-. 30373334E 

-02  • YD  EO . . 394 OE- 

04 

'IT. 

Cl HE  <SEC. > S 

PEED  <K  ' 3.' 

HEhDIH‘5<P. 

3i_T.  <KM.  < 

'3P . P!49  <KM  ■ 

1 

. 3000E+01 

. 1 349E-0 1 

. 1571E+01 

. 2 0 04E-0 1 

, 0 0 0 OE  + 0 0 

2 

. 3000E+0 1 

. 2774E-0 1 

. 1571E+0 1 

. 3 1 72E-01 

• OOOOE+OO 

3 

. 9000E+01 

.4279E-01 

. 1571S+01 

. 1 37 3E+00 

. OOOOE+OO 

4 

. 1200E+02 

. 5339S-01 

. 1533E+01 

. 3S93S400 

. 2277E-03 

5 

. 1500E+02 

. 7550E-01 

. 1542E+01 

. 5403E+00 

. 3423E-02 

6 

. 1300S+02 

. 9323E-01 

. 1503E+01 

. 7929E+00 

. 1543E-01 

. 2 1 00E+02 

. 1 122S+00 

. 1454S+01 

. 1 099E+01 

. 43S4E-01  ■' 

.3 

. 2400E+02 

. 1323E+00 

. 1393S+01 

. 1432E+01 

. 97  I'-E-h  1 /' 

3 

. 2700E+02 

. 1540 5*00 

. 1 333E+01 

. 1332E+01 

. 1343E+00 

10 

. 3000E+02 

.17735*00 

. 1275E+01 

.2331S+01 

. 3 1 5 1 E + o 0 

1 1 

• 3300E+02 

. 2024E+00 

. 1213E+0C 

. 2900E+0 1 

. 4939E+00  It 

12 

. 330 OS +02 

.2294E+00 

. 1 133E+01 

. 3501E+01 

. 7333E+00  . 

13 

. 3900E+02 

. 2533S+00 

. 1 1 15E  + 01 

. 4135E+01 

. 1 042E+01 

14 

. 4200S+02 

. 2390E+00 

. 1 033E+01 

. 4392E+01 

. 1 421E+01  *'■ 

1? 

. 4500S+02 

. 3209S+00 

. 1 0 1 7E+0 1 

. 5391 £+01 

. 1 3S3E+0 1 [v 

16 

. 4300E+02 

. 3544E+00 

. 9399E+00 

. 3530E+01 

. 2435E+01  r 

i? 

. 5 1 OOE+02 

. 390 IE* 00 

• '??3  '?E  ^ 0 0 

. 7435E+01 

. 3037E+01 

13 

. 54 00E+02 

. 4235E+00 

. 3793E+00 

. S397E+01 

. 3347E+0 1 V 

13 

. 57 00E+02 

. 4593E+00 

. 5370E+00 

. 9413E+01 

. 4727E+01  p 

20 

. 3000E+02 

. 5135S+00 

. 793  3£+0 0 

. 1 0495*03 

. 5735E+01 

21 

. 3300S+02 

. 5305E+00 

. 7573E+00 

. 1 132E+02 

. 3331  £ + 01  f ii 

. 33 00E+02 

. 31 07S+00 

• 72  0?E^00 

. 1230E+02 

.S177E+01  | :] 

23 

.3900S+02 

. 334 IE+00 

• 5550E+00 

. 1 4043+02 

. 9332E+01  H 

24 

• 72 00E+02 

. 7203E+00 

• 555  0E+00 

. 1 5 32E+02 

• 1 1 2^E+02 

25 

. 7500S+02 

. 7309E+00 

. 3220E+00 

• 1 666E+ 02 

. 1 303E+02 

23 

. 7S00E+02 

. 3444S+00 

. 5927E+00 

. 1 3 05E+ 02 

. 1503E+02  ' f 

27 

.31 00E+02 

.91 15E+00 

. 5352E+00 

. 1949E+02 

. 1723E+02  !l 

23 

. 34 00E+02 

. 9324E+00 

. 5393E+00 

. 2 093E ♦ 02 

. 1937E+02 

29 

.370 OS +02 

. 1057E+01 

. 5150E+00 

. 2252E+02 

• 22  ? OE  + Oc  f 

30 

. 9000E+02 

. 1 1 35E+01 

. 4921E+00 

. 24 l 1E+02 

. 2517E+02 

31 

. 9300E+02 

. 121  3E ♦ 0 1 

. 47  06E+ 00 

• 2574E  02 

. 2329E+02  fa 

32 

. 93 00E+02 

. 1305E+01 

• 45 04E  + 0 0 

. 2742E+02 

. 3133E+02 

33 

. 99 00E+02 

. 1393E+01 

. 43 1 4E+00 

. 29 1 5E ♦ 02 

. 35 31 £+02 

34 

. 1 020E+03 

. 1492E+01 

. 4133S+00 

. 3092E+02 

. 3924E+02 

35 

. 1 05 OS* 03 

. 1593E+01 

. 3274E+02 

. 4 34  7E+02 

33 

. 1 05 OS +03 

. 1345E+01 

. 331 3E+00 

. 3432E+02 

• 4'r’.  0 1 E ^ c. 

37 

.111 0E+03 

. 131 0E+01 

• 5555F+O0 

. 3354E+02 

. 5233E+02  ^ 

33 

. 1 140E+03 

. 1927E+01 

. 3529E+00 

. 3351 £+02 

• 55 1 OE + 02 

3? 

. 1 1 70S *03 

. 2050E+01 

• 34 OOF +00 

. 4053E+02 

• 6 36?E+0c! 

40 

..1200S+03 

. 2131E+01 

. 3230E+00 

. 4231E+02 

.6*6m<E  + 02  | 

IT. 

TIHE <3 EC.  • S 

PEED<K-S> 

HEhDI  H>5  • p> 

hLT. <KM.> 

■5P . PM6  ( K M • 1 

1 

.121 0E+03 

. 21 77E+01 

. 324 15*00 

. 4 3 31E  + 02 

. 7133E+02 

2 

. 1220E+03 

•2174E+01 

. 3202E+0O 

. 440OE+02 

. 7373E+02 

3 

. 1250E+03 

.2171E+01 

• 44»5:?E4‘  02 

• 75  75E^  ••  2 

4 

. 1240E+03 

. 2133E+01 

. 3 1 24E+ 0 o 

• 45  55E^  0.2 

. 7733E+02 

*, 

. 1250E+03 

. 2 1 35E  + 0 1 

• 3035E+00 

. 4301E+02 

• 755‘7E^  02 

if. 

. 1230E+03 

.2132E+01 

. 3045E+00 

• 4*£>i?f?E*0c* 

. 3 1 92E+02 

? 

. 1 27  OE ♦ 03 

. 2159E+01 

• 50  »*•*:•£  + 00 

• 4 7 5 i E + 02 

• 5 ?*->7E^  «.*2 

3 

. 1 23  OE  ♦ 0 3 

. 2 1 53E  + 0 1 

. w'H’sE^Om 

. 4794E+02 

. 3302E+02  < 

*3» 

. 1290S+03 

. 2153E+01 

. 2927E+00 

• 4c!57E^  02 

• 35  0*E + 02 

10 

. 1 300S+03 

.215  OE + 0 1 

. 2337E+0O 

.491 9E+02 

.401 1E+02 

IT. 

TIME  <?EC . > Z 

PEED  'K 

HE3DIH5  <P • 

ML T.  - HI.  • 

3P.Pl1i3’M1> 

1 

• 1 450E+03 

. 2 394E  ♦ 0 1 

. 2322E+00 

. 5792E+02 

. 1227E  + 0 3 

Wis  PAOi  is  USX  QUALITT  PRACTJCAM1 
J*0*  COrY  HNHSKSD  TO  ODQ  -60- 


. l*300E*03 

.2-37?E*01 

. 1322E*00 

■ 3?i3  9£*02 

. 1?9?E*03 

3 

. 17?0£*03 

. 2997E*01 

. 1 33?S*00 

. 7243E*02 

.301 OE ♦ 0 3 

4 

. 1900E+03 

. SSAAE^Ol 

.101 0E*00 

. 7S09E+02 

, 2473E*03 

? 

. 20S0£*03 

. 3S03£*0l 

. 33 3 IE- 01 

. 32A1E*02 

. 3007E*03 

(3 

. 2200E*03 

. 4 3 1 *3E  ♦ 0 1 

. 4337E-01 

. 3?99E*02 

. 3-30?£*03 

7 

. 23?0E*03 

. 4939S*0l 

. £3 0-3E-0 1 

. 3924E+02 

.423*E*03 

•? 

. 2?00E*03 

. ?72?E*01 

. 9443E-02 

. 5944S+03 

. ?074S*0  3 

9 

.2*?0£*03 

. -3737E*01 

1 304E-03 

. 39?7E*02 

. 599?£*03 

10 

. 290 OS* 03 

. 3424E+0 1 

1 7 12E-02 

.3A?9£*02 

. 7 1 09E  ♦ 03 

IMT. 

TIME* SEC. > 

SPEED**  S' 

HEAD IMA  *P • 

alt. *KM. > 

SP . PMS • * M • 

1 

. 2S10£*03 

. '34£4E>0 1 

-. l?41E-0£ 

. 3 9?3E ♦ o2 

. 7I9?E*03 

9 

. 2S20£*03 

. 3424Ef01 

137 OS- 03 

. 39?7E*03 

. 727?E*03 

3 

. 39 3 OS ♦03 

. S424S+01 

-. 1 IASS- 03 

. 39?-3£*02 

. 73?9S+03 

4 

.234 OS ♦03 

. 3434E+01 

1 037E-03 

. :39??E*02 

. 7442E*03 

? 

. 23?0E*03 

. 3424E*0l 

-. 3??9E-03 

. $9?4£*02 

. 7?2?E*03 

3 

. 23-3  OE*  03 

.3424S*0l 

■334AE-03 

. 39? 3£* 02 

. 7-303E*03 

. 2970E*03 

•3424S*0l 

-.?1 34E-03 

. 39? ?£*0£ 

. 7-39  IE* 0 3 

3 

. 3390E*03 

. 3434E+01 

-. 3433E-03 

. 39?3E*02 

. 7774E+03 

9 

. 3990E*03 

. 3424E+01 

171  OS- 03 

. 39?2S  ♦ 03 

. 73?7S*03 

10 
EMTEP 
- . 0 03 

. 2900E+03 
VALUES  FDP 

. 3424S*01 

CA.TL2. TLD 

.331 4E-0-3 

. 39?3E*02 

. 794 OH* 03 

1 0. 
jr( 

EMTEP  VALUES  FOP  MT.  IOP* ICE* IPP» IBS* JCV» IMA. JMS 
£♦013 

IMT.  TIMS*  SEC.  ••  SPEED  <!*'•$.>  HSADIMS*R*  ALT.  > * M.  • SP. PMS *K!1  > 

IMT.  TIME* SEC.  * -SPEED «.K  S'  HEAD I MS *°  • ALT.  *'KM.  ‘5P . P!1G * f* !1  • 

10  . 1 30O£*03  .30  3 3E*01  .?209E*00  .*3493£*02  . 7492E*o2 

EMTEP  VALUES  FOP  CA.  TL2*  TLD 
01 
13?. 

140. 

EMTEP  VALUES  FOP  MT • I OP.  ICE*  IPP*  IBS*  JCV.  IMA.  JHS 
£4212 

IMT.  TIME*  SEC.  SPEED**'  S*  HEAD  I MS  < P • ALT.  <s-M.  * 5P.PMS<KM* 

IMT.  TIME*  SEC.  » SPEED**  S*  HEAD  IMA  *P  • ALT.  **M.  * SP.PH9*TM» 

10  . 3900£*03  .30??E*01  -.9422E-01  . 114IE  + 03  .7139E*03 
EMTEP  VAL'JES  FOP  CA*TL2*TLD 
-.  003 
13?. 

140. 

EMTEP  VALUES  FOP  MT*  IQP*  ICE«  IPP*  ID'S.  JCV.  IMA.  JMS 
£ ♦£££ 


EMTEP  VALUES  FOP  PSSS*EPS 

0.  M 

. 000003 


10 

. 3900E*03 

.3  0-3  3£  ♦ 0 1 - . 

30-33E-01 

. 1 430E*0  3 

. 719'aE*03 

\M 

£0. 

-.  70333?£4E-03  . VO  E-?. 

-.  30-3  3E- 

01 

10 

. 3900E*03 

. 30-3  3E*  01  -. 

7 005E-04 

. 1??SE*0? 

. 7 1 9?E ♦ 0 3 

*M 

ES>. 

-.  7034-344 3E-03  . VO  ES. 

-. 7005E- 

04 

10 

. 3900S*03 

. 30»3£*0l  -. 

437 IE- 07 

. l??-3£*03 

. 719?E*03 

;m 

SO. 

-.  70343-437 

E-03  • VO  E9. 

-. 4371E- 

07 

IMT. 

TIME  <SEC . .' 

SPEED  <K> S > HEAD IMA  *P • 

ALT.  **  M.  ,• 

•5P.  PMS  <*<M  • 

1 

. 1 4?0£*03 

. £347E*01  . 

43 3 IE* 00 

. 3013E*03 

1029E+0? 

3 

. IS00E+03 

. 249'3E*0t 

399£E*0O 

. 94A^h*03 

. 1 347E+03 

3 

. 1 7?0E*03 

. 373*3E*01 

331 7E*00 

. 1 09 IE* 03 

. 171 1E*03 

4 

. 1900E*03 

. 3134E+01  . 

3*3  3?E  ♦ 00 

. 1230E+03 

.£13*3E*03 

5 

. 3 0?  OE  *03 

. 3?34£*01  . 

1973E*00 

. 1 333E+0? 

• 3-301E+0S 

£ 

. 330 OS *03 

. 4007F+0 1 . 

1 3??E*00 

. 1 423£*0  3 

. 314?E*03 

♦ 

. £3?O£*03 

. 4-304E*0 1 . 

3O39E-01 

. 1 49?E*03 

. 377  IE* 03 

■3 

. £?00E*03 

. ?373E ♦ 0 1 

3333E  — 0 1 

. 1?37E*03 

. 4493E*03 

■ ? 

. 23? OS +03 

.3427E+01  . 7 1 17E-02 

. 1???E>03 

. ? 3?3E  + 03 

10 

. 2300E+03 

. 3033S+01  -. 7397S-03 

. l??3E«-03 

. 3407E+03 

IMT. 

TIME ( SEC. » 

SPEED  (K/S.  HEaDIMG (P • 

atT. (KM. ) 

GP . PMG  < KM> 

1 

. 231 0E+03 

. 3033E+01  -. 3333E-03 

. 1 3G3E+03 

. 3433E+03 

2 

. 2320E+03 

.3033E+01  -. 3073E-03 

. 1??3E*03 

. 3?3?E>03 

3 

. 233 OE* OS' 

. 3033S+01  -. 5313E-03 

. 13?3E«-03 

. 3344E+03 

4 

. 2340E+03 

.3033S+01  -.4??3£-03 

. l??3E+03 

.3722E^03 

? 

•2330E+03 

. 3 03  3E*  0 1 - . 3799E-0 3 

. l??3E+0? 

. 3301E+03 

3 

.233 OS +03 

.3033E+01  -. 3039E-03 

. 1 ?53E^03 

. 3330E^03 

.2370S+03 

•3033E+01  -. 2279E-03 

. 1??3E*03 

. G953E+03 

3 

. 2330E+03 

. 3033E+01  -. l?19E-03 

. l??3£+03 

. 7037E+03 

9 

. 2390E+03 

. 3033E+01  -. 7?97E-04 

. l??3E+03 

. 7 1 1GE+03 

10 

EMTEP 

-.OO? 

. 2900E+03 
VPL'JES  FQP 

. 3033E+0 l . 229?E-03 
CP. TLZ. TLD 

. 1??3E>03 

. 7 19?E+03 

to. 
a?. 

EMTEP 
2401.2 
IMT. 
IMT. 

10 
EHTEP 
01 
13?. 
140. 
EMTEP 
24212 
IMT. 
IHT. 

10 
EMTEP 
013 
13?. 

140. 
EMTEP 
24222 
EMTEP 
0. 

. 000002 
10 

XM  SO. 

10 


VALUES  POP  MT. IQP. ICE. IPP. IBG. JCV. IMG. JMS 


TIM£<SEC.  • 
TIME (SEC.  ' 

. 1300E+O3 

vaL'JES  CQP 


SPEED  <KV'S> 
SPEED  < S j 
. 1973E*01 
CP. TLZ. TLD 


HE9DIMG  (P> 
HEMDIMG • P 1 
. 3??9E*00 


PLT. <KM. > 
3L.T.  (KM.  > 
. 7227E+02 


VSL'JES  FQP  MT. IQP. ICE. IPP. IBG. JCV. IMa* JMS 


TIME  (SEC.  > 
TIME (SEC.  • 
. 2900E+03 
VGL'JSS  FCP 


SPEED  (KV  3 • 
SPEED (K  S • 
. 73 OSS ♦01 

CP. TLZ. TLD 


ME3DIMG  (P> 
HEaDIMG'P' 
. 27 14E-01 


vauJES  POP  MT.  IQP.  ICE. IPP. IBG. JCV. IMP. JHS 


aL'JES  POP 


2900E+03  . 7773E+01  -.7393E-0I 

1030?297S-01  . VQ  E9.  -.7393E- 
2'? 0 OE ♦ 03  . 7301E+01  . ?990E—  04 


KM  E9.  -.  10S07070E-01  . VQ  E 


IMT. 

I 

3 

3 

4 
? 
3 


3 

4 
10 

IMT . 
1 
3 

3 

4 
? 

£ 


TIME ( SEC. > 
. 14? OS* 03 
. 1300E+03 
. 17? OS +03 
. 1900E+03 
.20?0E+03 
. 2200E+03 
.23?0E+03 
.2?00S+03 
.23?0£+03 
. 2300E+03 
T I ME  < SEC..' 
.2310E+03 
.2320S+03 
. 23  3 0E+ 0 3 
. 2340E+03 
. 23?0E+03 
. 2330E+03 


SPEED  (K/S> 
. 2170E+01 
. 2400E+01 
. 2337E+01 
. 2973S+01 
. 33?0E*0l 
. 3303E+01 
.4373E+01 
.?123S+01 
. 3 137E+01 
. 7301E+01 
SPEED  <*/ S • 
. 7301E+01 
. 7301S+01 
. 7301E+01 
. 7S01S+01 
. 7301E+01 
. 730 1E+0 1 


9.  . ?930E- 

MEar i mg  < p • 

. 399 1 S+ 00 
. 5233S+00 
. 44S1E+00 
. 3324E+00 
.2737E+00 
. 1924E+00 
. 1 139S+00 
. ?437E-0l 
. 1 139E-01 
31 33S-04 
MEPDIMG • P ' 
-.232  IE- 04 
-. 2303S-04 
-. 2 13  IE- 04 
-. 1373E-04 
-. 1331E-04 
-. 1243S-04 


GP.  PMG  (KM  < 
GP . PMG 'KM' 
. 3???E+02 


3LT. <KM. > 
aLT.  (KM.  • 
. 2000E+03 


GP . PMG CKM . 
GP . PMG (KM • 
. 3393E+03 


. 1374E+03 
01 

. 1333S+03 
04 

PlT.  (KM.  ,■ 
. 9030S+02 
. 1 033E+03 
. 1233S+03 
. 1432E+03 
. 1331E+03 
. 1704E+03 
. 1797E+03 
. 1333E+03 
. 1332E+03 
. 1SS3S+03 
9lT.  (XM.  • 
. 1333E+03 
. 1 333E+0  3 
. 1 333E+0  3 
. 'I3G3E+0? 
. 1 333E+03 
. I 333E+03 


. 3390E+03 
. 37  01E  + 03 

GP . PMG (KM • 
. 90  34E+02 
. 1 133E+03 
. 1?1 7S  + 03 
. 1393E+03 
.2337S+03 
. 2S43E+03 
. 34331+03 
. 4120E+03 
. 4937E+03 
. 3944E+03 

GP . PMG ( K M • 
. 3019E+03 
. 3093E+03 
•3171E+03 
. 32  47S+03 
. 332  3S+0.3 
. 3 393E+03 


-62- 


. 2S70E+03 

. 73  0 1 E ♦ 0 1 

-.931.2E-05 

. 1995E+03 

. 6474E+0? 

3 

.iS'SOE  + OS 

. 7S01E+01 

■sl!i5£“0? 

. 189“5F>|*»? 

. •i5,5  0E*  03 

3 

.3S*0E+03 

. 7301E+01 

-.301 3E-05 

• 1 995E+ 03 

. 6626E+03 

10 

•2900E+03 

. 7301E+01 

. 1 3* IE- 06 

• 1 395E ♦ 0 3 

. 67  01E  + 03 

ATS9 

VSL'JES  £□<? 

Cfl«  T!_2»  TLB 

• K • F 


( 


U 

' 


DEPARTMENT  Of  THE  NAVY 


